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A COMPUTER PROGRAM FOR PLOTTING STRESS -STRAIN DATA FROM 
COMPRESSION, TENSION, AND TORSION TESTS OF MATERIALS 

By Anne Greenbaum and Donald J. Baker 
Langley Directorate, U.S. Army Air Mobility R&D Laboratory 

John G. Davis, Jr. 

Langley Research Center 

\ 

SUMMARY 

A computer program for plotting stress -strain curves obtained from compression 
and tension tests on rectangular (flat) specimens and circular -cross -section specimens 
(rods and tubes) and both stress -strain and torque -twist curves obtained from torsion 
tests on tubes is presented in detail. The program has the capability of plotting individ- 
ual strain -gage outputs and/or the average output of several strain gages and the capabil- 
ity of computing the slope of a straight line which provides a least -squares fit to a speci- 
fied section of the plotted curve. In addition, the program can compute the slope of the 
stress -strain curve at any point along the curve. The program, its subroutines, and their 
variables are listed and defined. The computer input and output for three sample prob- 
lems are presented in printed and plotted form to aid the user in setting up and utilizing 
the program. 

The program is written in FORTRAN IV language for the Control Data 6000 series 
digital computer with the SCOPE 3.0 operating system and requires approximately 
110000 octal locations of core storage. A typical problem with seven curves to be plotted, 
each curve containing approximately 200 data points, requires approximately 50 seconds 
of central processing time. Output from the program is stored on tape and then used by 
CalComp or Gerber plotters to construct curves. The input data may be stored on mag- 
netic tape, a data cell, or punched cards. 

INTRODUCTION 

Reduction of test data to a presentable form, usually hand -drawn plots, can be a 
time-consuming task. This is especially true if a large number of data points are to be 
plotted or if the data must be combined in some manner or multiplied by scaling factors 
prior to plotting. A specific example is construction of the shear -stress— strain curve 
from strain-gage data obtained from a torsion test on a tubular specimen. Since there 



has been a continuing requirement at the Langley Research Center to make such data 
reductions, a computer program was developed to automate the task of constructing 
stress -strain plots for compression and tension tests on rectangular (flat) specimens and 
circular -cross -section specimens (rods and tubes) and stress -strain and torque -twist 
plots for torsion tests on tubular specimens. The program plots individual strain-gage 
outputs and/or the average output of several strain gages. The program also has the 
capability of computing the slope of a straight line which provides a least-squares fit to a 
specified section of the plotted curve and can compute the slope of the stress -strain curve 
at any point along the curve. The program is written in FORTRAN IV language for the 
Control Data 6000 series digital computer with the SCOPE 3.0 operating system. The 
input data may be stored on magnetic tape, a data cell, or punched cards. The equa- 
tions have been programed so that either the International System of Units or the U.S. 
Customary Units may be used. Conversion factors relating the two systems are given in 
reference 1, and those pertinent to the present paper are presented in appendix A. This 
paper describes the program and presents sample problems to aid the user in setting up 
and utilizing the program. 


SYMBOLS 

£ length, meters (inches) 

P load, newtons (pounds) 

R radius, meters (inches) 

T torque, newton -meters (inch-pounds) 

t thickness, meters (inches) 

w width, meters (inches) 

5 displacement, meters (inches) 

e strain 

9 twist, radians (degrees) 

o’ stress, pascals (pounds force per square inch) 
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t shear stress, pascals (pounds force per square inch) 

Abbreviations: 

KSI kilopounds force per square inch 

LVDT linear variable differential transformer 

MPA megapascals 

MICRO IN/IN micro inches per inch (see description of XMULT in appendix C) 
MICRO M/M micro meters per meter (see description of XMULT in appendix C) 
O.D. outer diameter, meters (inches) 

PA pascals (newtons per square meter) 

PSI pounds force per square inch 

SG strain gage 


PROBLEM DESCRIPTION 

The task described herein is the development of a computer program to construct 
plots of stress -strain data from compression and tension tests and to construct plots of 
stress -strain and torque -twist data from torsion tests. At the beginning of the effort, th« 
following criteria on program capability were selected: 

(1) The program must be capable of reading load, strain, and displacement input 
data, which are arranged sequentially with respect to time and expressed in appropriate 
engineering units, from magnetic tape, a data cell, or punched cards. 

(2) The program must be capable of plotting data in either the International System 
of Units or U.S. Customary Units. Data recorded in either system of units can only be 
plotted in the respective units. 

(3) The program must be capable of constructing compression or tension stress - 
strain plots for either rectangular (flat) specimens or circular-cross-section specimens 
(rods and tubes) and constructing shear -stress— strain and torque -twist plots for thin- 
walled circular tubes. 
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(4) For compression or tension tests, the program must be capable of constructing 
plots of stress as a function of the output of individual strain gages and/or the average 
output of up to four specified strain gages (or other strain transducers). 

(5) For torsion tests, the program must be capable of constructing plots of stress 
as a function of the output of individual strain gages, as a function of the average absolute 
output of the ±45° strain gages in a strain rosette, or as a function of the average absolute 
output of all ±45° strain gages in up to six strain rosettes. Also, the program must be 
capable of constructing torque -twist plots by utilizing the output of LVDT’s or other 
transducers which measure displacements associated with the rotations at two locations 
along the length of a tubular specimen. 

(6) The program must be capable of computing the slope of a straight line which 
provides a least-squares fit to a specified section of each averaged curve (stress -strain 
or torque -twist) and computing the tangent modulus (slope) at all points along the average 
stress -strain curve. 

(7) The program must be capable of omitting from the calculations and plots any 
gages which the user specifies, and it must be capable of making more than one type of 
plot from one input of the same data set. 

These criteria have been met and table I summarizes the plotting cases that have 
been programed. Ten cases have been programed and the table indicates specimen and 
test type, the quantities plotted on each axis, and whether moduli or slope calculations 
can be made for the case. Figure 1 shows sketches of specimen geometry for compres- 
sion, tension, and torsion test specimens. 

PROGRAM DESCRIPTION 

The computer program is written in FORTRAN IV language for the Control Data 
6000 series digital computer with the SCOPE 3.0 operating system. With the present 
dimensions, the program requires approximately 110000 octal locations of core storage. 

A typical case with seven curves to be plotted, each curve containing approximately 
200 data points, requires approximately 50 seconds of central processing time. During 
job execution, input data read from a data cell or magnetic tape are temporarily stored 
on a disk file (designated as tape 30 in the program listing). As currently written, the 
program is limited to a maximum of 20 data channels (18 if stored on a data cell) and 
414 data frames per run, with 398 frames per run plotted. Dimension statements can be 
changed to allow for a larger number of channels and frames. A data channel is a listing 
that contains all data values derived from one output device (force transducer, strain gage, 
LVDT, etc.). A data frame is defined as all data values recorded at one value of time. 
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Thus, one data frame contains a data value for each output device. To use the program 
as it is currently written, the data must be identified by frame number and stored sequen- 
tially by frame number. 

The main program, identified herein as BECKPLT, and its subprograms are iden- 
tified and discussed. Except where noted, all subprograms were written for the task 
described herein. Appendixes B, C, and D give a source program listing, definitions of 
the FORTRAN variables, and the Langley Library subroutine MATIN V and listing, respec- 
tively. A flow chart is shown on pages 6 to 8. The flow chart indicates that, first, 
instructions that indicate the title of the test, the type of plotter to be used (Gerber or 
CalComp), how the data are stored (magnetic tape, punched cards, or data cell), the spec- 
imen geometry (flat or tube), and test type (tension, compression, or torsion) are read 
from cards 1 and 2. For simplicity, the type of specimen is given in the program as 
flat or tube (includes both rods and tubes). The Gerber plot is a smooth continuous 
curve, whereas the CalComp plot is incremented in small steps. The CalComp plotter 
requires less time to develop a curve than a Gerber plotter, but the curve is of lower 
quality. Then the instructions that indicate the particular experiment being analyzed 
(referred to hereinafter as run), the number of data channels and frames to be plotted, 
and the curves to be plotted (refer to table I) are read from cards 3 and 4. The data are 
read from the designated storage device (tape, cards, or data cell) and then the specimen 
geometry is read from instruction card 5. Next, the load data are converted to stress 
values by multiplying by the appropriate constants which are a function of test type and 
specimen geometry. These stress values are stored in the Y -array. One subroutine is 
next selected to compute and store the strains in the X-array. Selection is based on the 
type of test and specimen. Possible choices of subroutines are: COMPIND (compression 
or tension test, individual gage plots), COMPAVG (compression or tension test, average 
plot), TORIND (torsion test, individual gage plots), and TORAVG (torsion test, average 
plot). At this point, the coordinates for the stress -strain curve are computed and stored 
on tape for subsequent use by the selected plotter. Several subroutines can now be 
requested in the program. Subroutine SLP can be called to calculate the slope of a 
straight line which provides a least-squares fit to a specified portion of the stress-strain 
curve. If requested by user instructions to compute the tangent modulus (slope) of the 
stress -strain curve as a function of stress, subroutine PROPLIM is called. Subroutine 
PROPLIM computes the coefficients of a second -degree polynomial that provides a least - 
squares fit to the stress-strain curve. Next, the slope of the stress-strain curve is 
approximated by the first derivative of the polynomial. The user is cautioned that mean- 
ingful tangent-moduli results will be produced only when the second-degree polynomial 
is a good approximation to the stress -strain curve. Preliminary attempts to compute 
the tangent moduli by piecewise fitting the data (5 to 15 points) with a second-degree 
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polynomial generated plots (tangent moduli as a function of stress) of no practical value. 
Variation of the computed values of tangent moduli ranged up to 15 percent of the mean 
value, and the curves subsequently plotted had a sawtooth appearance. 

More than one plot (case) can be developed from a data set (run) by utilizing 
instruction card 6. This card indicates the frames to be used and the type of additional 
plots to be generated. An example of how this feature of the program might be used is 
the construction of the average shear -stress— strain and torque -twist curves (case 10 in 
table I) after the program has constructed the individual stress -strain curves (case 7). 

With one computer job submittal the stress -strain and/or torque -twist curves 
for each specimen of a group of specimens can be developed by repeating instruction 
cards 3 to 6 for each specimen. 

The program just described was written primarily to develop data plots for the 
cases listed in table I. However, this program can be modified so that it can be used to 
construct data plots for other types of engineering tests. These modifications would 
include the subprograms that label axes and transform data. The location of the changes 
can be determined by following the flow chart just discussed. 

PROGRAM USAGE 

A typical program deck setup, with punched card input, is shown in figure 2. Sys- 
tem control cards are utilized to position magnetic -tape or data -cell inputs for use by 
the source program. User instruction cards, described later in this section, are uti- 
lized for control of plots and computation of the slope of the curves. A detailed descrip- 
tion of the output is also given in this section. In order to keep user instructions at a 
minimum, a pattern for identifying the strain gages and rosettes was selected. (Refer to 
appendix E.) 

Input 

Data format. - Input data for the program may be recorded on magnetic tape, a data 
cell, or punched cards by frame. Each frame includes the load and all strain -gage data 
recorded at a particular time. Frames are identified by integers which increase mono- 
tonically with increasing time from the start of the test. Load or torque values must be 
stored in the first data channel. 

If a magnetic tape is used, test data for each frame must be stored on the tape in 
the following format which is currently used at Langley Research Center: 
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Word 

Description 

Type 

1 

Mode 

Integer 

2 

2 

Integer 

3 

3 

Integer 

4 

Frame 

Integer 

5 

5.0 

Float 

6 

6.0 

Float 

7 

7 

Integer 

8 

8 

Integer 

9 

9 

Integer 

10 

10 

Integer 

11 

11.0 

Float 

12 

12.0 

1 

Float 

13 

13.0 

Float 

14 

14.0 

Float 

15 

15.0 

Float 

16 

Load or torque 

Float 

17 

1 

Data channels 

Float 

t 

114 




Record length is 114 words. At the end of a run, a record is written with Mode = 8888. 
At the end of all data on a tape, a record is written with Mode = 9999. There is an EOF 
after this record. Mode = 2 for all other data frames. 


When a data cell is used, input for each data frame must be in the following format: 


Line 

FORTRAN variable name 

Description 

Format 

i 

ITST, MRN, MD, IFRA 

Test, run, mode, frame 

3112, 24X, 112 

2 

DTA(IFRA,1),DTA(IFRA,2), 
DTA(IFRA ,3) ,DTA(IFRA,4) , 
DTA(IFRA,5),DTA(IFRA,6) 

Load or torque, first 
five data channels 

6E12.4 

3 

DTA(IFRA,7),DTA(IFRA ,8) , 
DTA(IFRA,9) ,DTA(IFRA,10) , 
DTA(IFRA ,1 1) ,DTA(IFRA ,12) 

Next six data channels 

6E12.4 

4 

DTA(IFRA,13),DTA(IFRA,14), 

DTA(IFRA,15),DTA(IFRA,16), 

DTA{IFRA,17),DTA(IFRA,18) 

Next six data channels 

6E12.4 
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(For a test with more than 18 channels stored on a data cell, the READ statements in 
subroutine RDCELL must be changed.) 

Input data on punched cards must be in the following format for each data frame. 


Card 

FORTRAN variable name 

Description 

Format 

1 

IFRA, DTA(IFRA,1), 
DTA(IFRA ,2) ,DTA(IFRA,3) , 
DTA(IFRA ,4) ,DTA(IFRA ,5) , 
DTA(IFRA,6),MD 

Frame, load or torque, 
first five data 
channels, last card 

14, 6E12.4, 14 

2 

DTA(IFRA,7) ,DTA(IFRA,8) , 
DTA(IFRA,9) ,DTA(IFRA,10) , 
DTA(IFRA,11),DTA(IFRA,12), 
MD 

Next six data channels, 
last card 

6E12.4, 14 


Repeat card 2 for all channels. For all data cards except the last card, MD = 2. On last 
card, MD = 8888 (end of data). 

User instruction cards .- The following cards are supplied by the program user to 
control the graphic and printed output. Cards 1 to 3 have the same format for tension, 
compression, or torsion tests. However, the format for cards 4 to 6 is a function of test 
type as indicated herein. 

User instruction cards 1 to 3 follow: 


Card 1 Test title (4A10) 

Column 

1 to 40 Title of test. Any characters in columns 1 to 40. This title is printed 

before all other output and is plotted in front of the graphs. (See 
section entitled "Sample Problems.") 


Card 2 Plotting control (914) 


Column FORTRAN variable Value 

4 IPLT 1 

0 


Use Gerber plotter 
Use CalComp plotter 


11 




Column 

8 


FORTRAN variable 
IPEN 


Value 


12 

16 

20 

21 to 24 


25 to 28 
29 to 32 
36 


INPT 

ISPEC 

ITEST 

IXSIZE 


IYSIZE 


IOMOVE 


IUNITS 


1 Use Leroy pen (for CalComp 

plotter only) 

0 Use ballpoint pen (for CalComp 

plotter only) 

0 For Gerber 

0 Input data on magnetic tape 

1 Input data on punched cards 

2 Input data on data cell 

0 Flat specimen 

1 Tube specimen 

0 Tension or compression test 

1 Torsion test 

Length of X-axis for each plot 
(in inches). For tension or 
compression tests, IXSIZE 
= (IOMOVE) x (number of curves 
+ 1). For torsion tests, IXSIZE 
= (IOMOVE) x (number of curves 
+ 1) + (number of rosettes). The 
last term accounts for the extra 
space between plots for different 
rosettes 

Length of Y-axis in inches. If this 
field is blank or = 0, length of 
Y-axis defaults to 9 inches 

Distance between plots. If this 
field is blank or = 0, defaults to 
2 inches 

3 Input and output data expressed in 

U.S. Customary Units 

l Input and output data expressed in 

SI Units 
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Card 3 Run title 


(4A10) 


Column 

1 to 40 Title of run. Any characters in colums 1 to 40. This title is printed 

before other output for this run and is plotted in the lower left-hand 
corner of the graph for this run. (See section entitled ’’Sample 
Problems.”) 

User instruction cards 4 to 6 for tension or compression tests follow: 


Card 4 

Run information 

(413 ,412, 313, 1812, 14X, 11) 

Column 

FORTRAN variable Value 


2 to 3 

JAC 


Number of channels (=1 + number 




of strain channels) 

4 to 6 

NRL 


First frame to be plotted 

7 to 9 

NRU 


Last frame to be plotted 

12 

INDAVG 

• 

Type of plot (see table I) 



0 

Case 1 or 4 



1 

Case 2 or 5 



2 

Case 3 or 6 

13 to 20 

IGAGES(l) 


Channels to be averaged. See 


IGAGES(2) 


sample problem 2 for example of 
input data 


IGAGESt3) 




IGAGES(4) 



23 

MOD 

0 

Do not calculate slope 



1 

Calculate slope 



2 

Calculate tangent moduli 

24 to 26 

NRLSLP 


First frame to be used in calculat 


ing slope. May be omitted if 
MOD * 1. If MOD = 1 and this 
field is blank or 0, defaults to 
first frame plotted 
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Column FORTRAN variable Value 

27 to 29 NRUSLP 


30 to 65 NOPLT(l) 

NOPLT(2) 


Last frame to be used in calculat- 
ing slope. May be omitted if 
MOD * 1. If MOD = 1 and this 
field is blank or 0, defaults to 
last frame plotted 

Strain gages to be omitted from 
plot 


66 to 79 
80 


NOPLT(18) 

MORE 


0 No additional plots (cases) 

requested for this data set 
(card 6 must be omitted) 

1 Additional plots to be made for 

this run 


If data are on punched cards, place data for this run after card 4. 


Card 5 Sectional data (4F12.0) 


Column 

FORTRAN variable 

Value 


1 to 12 

THICK 


Thickness of flat specimen, wall 
thickness of tube or radius of 
rod (circular cross section) 

13 to 24 

WIDTH 


Width of flat specimen; or outer 
diameter of tube specimen or 
rod 

25 to 36 

SCFAC 

+1 

-1 

Loads are positive 

Loads are negative 

Loads are multiplied by this 
variable before calculating 
stresses 
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Column 


FORTRAN variable 


Value 


37 to 48 RMAXSN 


Maximum strain value to be 
plotted 


If additional plots were requested for this run, place the following card or cards after 


card 5. 




Card 6 

New set of run information 

(313, 412, 313, 1812, 17X,I1) 

Column 

FORTRAN variable 

Value 


1 to 3 

NRL 


First frame to be plotted 

4 to 6 

NRU 


Last frame to be plotted 

9 

INDAVG 


Case type (see table I) 



0 

Case 1 or 4 



1 

Case 2 or 5 



2 

Case 3 or 6 

10 to 17 

IGAGES(l) 

IGAGES(2) 

IGAGES(3) 

IGAGES(4) 


Channels to be averaged 

20 

MOD 

0 

Do not calculate slope 



1 

Calculate slope 



2 

Calculate tangent moduli 

21 to 23 

NRLSLP 


First frame to be used in calculat- 
ing slope. May be omitted if 
MOD * 1. If MOD = 1 and this 
field is blank or 0, defaults to 
first frame plotted 

24 to 26 

NRUSLP 


Last frame to be used in calculat- 


ing slope. May be omitted if 
MOD * 1. If MOD = 1 and this 
field is blank or O', defaults to 
last frame plotted 
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Column FORTRAN variable Value 

27 to 62 NOPLT(l) 

NOPLT(2) 


Strain gages to be omitted from 
plot 


63 to 79 
80 


NOPLT(18) 

MORE 


0 No additional plots (cases) 

requested for this data set 
(card 6 must be omitted) 

1 Additional plots requested for this 

run 


User instruction cards 4 to 6 for a torsion test follow: 
Card 4 Run information (8I3,18I2,19X,I1) 


Column 
2 to 3 


4 to 6 
7 to 9 
12 
15 


18 


FORTRAN variable 
JAC 


NRL 

NRU 

NROS 

INDAVG 


MOD 


Value 


0 

1 

2 

3 

0 

1 

2 


Number of channels (= 1 + number 
of strain channels + number of 
deflection (6) channels (see 
fig- 1)) 

First frame to be plotted 
Last frame to be plotted 
Number of rosettes 
Case type (see table I) 

Case 7 
Case 9 
Case 10 
Case 8 

Do not calculate slope 
Calculate slope 
Calculate tangent moduli 
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Column 

FORTRAN variable 

Value 


19 to 21 

NRLSLP 


First frame to be used in calculat- 
ing slope. May be omitted if 
MOD * 1 . If MOD = 1 and this 
field is blank or 0, defaults to 
first frame plotted 

22 to 24 

NRUSLP 


Last frame to be used in calculat- 
ing slope. May be omitted if 
MOD * 1. If MOD = 1 and this 
field is blank or 0, defaults to 
last frame plotted 

25 to 60 
61 to 79 

NOPLT(l) 

NOPLT(2) 

NOPLT(IS) 


Strain gages to be omitted from 
plotting and averaging 

80 

MORE 

0 

1 

No additional plots (cases) 
requested for this data set 
(omit card 6) 

Additional plots requested for this 
run 

If input data are on punched cards, place data for this run after card 4. 
Card 5 Sectional data (6F12.4) 

Column 

FORTRAN variable 

Value 


1 to 12 

THICK 


Thickness of tube 

13 to 24 

WIDTH 


Outer diameter of tube 

25 to 36 

RLENGTH 


Distance between arm stations 


(see fig. 1) 


Length of R (see fig. 1) 


37 to 48 


ALENGTH 


Column FORTRAN variable 

49 to 60 SCFAC 


61 to 72 RMAXSN 


Value 

+1 

Loads are positive 

-1 

Loads are negative 


Loads are multiplied by this 
variable before calculating 
stresses 

Maximum strain value to be 
plotted 


If additional plots were indicated for this run, place the following card or cards after 
card 5. 


Card 6 New set of run information (713, 1812, 22X, II) 


Column 
1 to 3 
4 to 6 
9 

12 


15 


16 to 18 


FORTRAN variable Value 

NRL 
NRU 
NROS 
INDAVG 

0 

1 

2 

3 

MOD 0 

1 

2 

NRLSLP 


First frame to be plotted 
Last frame to be plotted 
Number of rosettes 
Case type (see table I) 

Case 7 
Case 9 
Case 10 
Case 8 

Do not calculate slope 

Calculate slope 

Calculate tangent moduli 

First frame to be used in calculat 
ing slope. May be omitted if 
MOD * 1. If MOD = 1 and this 
field is blank or 0, defaults to 
first frame plotted 
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Column 


FORTRAN variable 
NRUSLP 


Value 


19 to 21 


22 to 57 NOPLT(l) 


Last frame to be used in calculat 
ing slope. May be omitted if 
MOD * 1 . If MOD = 1 and this 
field is blank or 0, defaults to 
last frame plotted 

Strain gages to be omitted from 
plot 


NOPLT(2) 


No additional plots (cases) 
requested for this data set 

Additional plots requested for this 
run 

Repeat card 6 for each plot after the first plot to be made for this run (until 
MORE = 0). Repeat cards 3 to 6 for each data set in a computer run. 

Output 

Printed. - Examples of printed and plotted output are presented in the discussion of 
sample problems. The test title supplied by the user is printed first, followed by the test 
and specimen type and the plotter. Next, the first run title supplied by the user is printed, 
followed by a list of data from the first three channels of every fifth or tenth frame (every 
tenth frame if the total number of frames is greater than 300). This data list is followed 
by the specimen thickness and width (for flat specimen) or outer diameter (for tube speci- 
men). Next the coordinates of each stress -strain or torque -twist curve to be plotted are 
printed. These points are printed in pairs with the Y -value (stress or torque) first and 
the X-value (strain or twist) second. If tangent moduli are computed, they are listed 
next along with their associated stress -strain points. At the end of the printed output, 
the number of graphs plotted and the size of the axes are listed. 

Plotted.- The test title is printed vertically on the plotting paper before the first 
plot is~constructed. (See fig. 3.) The run title is printed in the lower leftrhand corner 
of the graph, below the X-axis. In graphs with more than one plot, the X-axis is labeled 


58 to 79 
80 


NOPLT(18) 


MORE 
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with a drafting dimension: -Jscale factor|— . In graphs with only one plot, both the X- 
and Y-axes are numbered at each tick mark. A grid of +'s is also drawn on the graph. 

Above the graph, the gage number and the frames used for each plot are written. 

For plots where the gages to be averaged are specified by the user, the program prints 
AVG OF and lists these gages. For plots where the gages to be averaged are fixed 
(cases 9 and 10 in table I), it prints AVG PLOT. 

If the tangent moduli have been calculated, they will be plotted as a function of stress 
on a 9- by 9-in. graph. On the X-axis will be STRESS in PSI or PA and on the Y-axis will 
be TANGENT MODULI in PSI or PA. 


Diagnostics 

This program contains a series of error diagnostics as follows: 

Fatal errors - 

(1) INVALID TEST NUMBER, ITEST = (test number specified). Neither tension, 
compression, or torsion test is indicated. (ITEST is not 0 or 1.) 

(2) INVALID SPECIMEN NUMBER, ISPEC = (specimen number specified). Neither 
flat nor tube specimen is indicated. (ISPEC is not 0 or 1.) 

(3) CANNOT HANDLE TORSION TEST ON FLAT SPECIMEN 

(4) TOO MANY CHANNELS, JAC = (number of channels specified). The number of 
channels specified is greater than 20. 

Nonfatal errors - 

(1) LAST FRAME TO BE PLOTTED LESS THAN FIRST FRAME 
FIRST FRAME = (first frame) LAST FRAME = (last frame) 

WILL GO TO NEXT PLOT 

This can be used to suppress plotting for a particular run, if desired. 

(2) FIRST FRAME FOR SLOPE GREATER THAN LAST FRAME FOR SLOPE. 
FIRST FRAME = (first frame for slope) LAST FRAME = (last frame for slope) 
WILL NOT FIND SLOPE. 

The program will print and plot the points for this run but will not calculate the slope of 
the line of best fit to any section of the curve. 

(3) FIRST FRAME FOR SLOPE INVALID. 

FIRST FRAME = (first frame specified for slope) 

WILL USE FIRST FRAME PLOTTED. 
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(4) LAST FRAME FOR SLOPE INVALID. 

LAST FRAME = (last frame specified for slope) 

WILL USE LAST FRAME PLOTTED. 

SAMPLE PROBLEMS 
Problem 1 

This example illustrates the input, output, and typical stress-strain plots obtained 
from a compression test on a tube. The tube was instrumented with nine strain gages. 
The output of all nine gages as well as the average output of gages 2,5, and 8 is plotted 
in figure 3. Input and output data are expressed in SI Units. The CalComp plotter with 
a Leroy pen was utilized to plot figure 3. A listing of the case instruction cards follows: 


/ . oqiposi . 0766 

r -1 1 

. 0 ij 1 0 1 A 
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(Card 5) 
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(Card 1) 
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The data set was obtained from a test on a graphite -epoxy tube. 
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The output listing for problem 1 follows: 


PROBLEM 1 C OPPRESSION T f S T 

COPPPFSSI ON OP T[wi; jr;^ y !- S 1 )\i TUH t: FP6CM6N 

US T MG CALCOMP PinTT^ wit'l LK'-’OY PF ‘-J 

TURF N'ljyocp 2 


FPflMP L'lAO 

6 *.l76)r«-o? i. 

10 ?.??<:. 7f+o.R 4. 

15 4.19R3F+03 4. 

2 0 ft. O'* RTF + 03 1 . 

7 5 P . 1 1 n 0 F + OR 1. 


SO-1 S>’ 

6341 r-f 5 3. Pi 3526-05 

14376-05 1.10056-04 

R 3ft ?f-:) S 1.9 39RP-04 
17126-05 2. 65976-04 

R 0 0 ft f - 0 4 7.53MF-04 




1 55 
160 
165 
17 0 
175 

THICKNESS = 


l .66! 5F *04 
1 . 7 43 RE *04 
1 .03016+04 
7.1099F+04 
2.23346+04 
l .3091F-03 


R.’l 78 c -r>4 
?. 4679F-04 
2. 7 013=- 04 
2 .92666-04 
3.0632F-O4 
OUT FR 


6 . 73006-04 
7.5R04F-04 
B .30416-04 
9 .0545F-04 
9 .5 3R 16-04 
OIAHETER = 


7. 6600F-02 
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Problem 2 


This example illustrates the input and output data from a compression test on a 
graphite -epoxy tube. The data are expressed in U.S. Customary Units. This problem 
illustrates the option where more than one type of plot is obtained from a data set (run). 
First, the output of nine strain gages and the average of the output of gages 2,5, and 8 
are plotted in figure 4. The tangent modulus for the averaged curve (gages 2, 5, and 8) 
is plotted in figure 5 as a function of applied stress. The dashed curve in figure 5 is 
hand drawn and represents the slope of the straight line which provides the least-squares 
fit to the data. Figures 4 and 5 correspond to case 5 in table I. The additional output 
requested for this computer run is a separate plot of the average output of gages 2,5, 
and 8 in figure 6. The slope of the stress -strain curve shown in figure 6 is 10.6 x 10 6 , 
as indicated in the printout and by the dashed line in figure 5. A listing of the input data 
follows: 



I B ifti 0 1 0 0 0 0 0 0 FORTRAN STATEMENT 


\ (Card 6) 

i ^ (Card 5) 
Tjx^ (Card 4) 
(Card 3) 


A 


(Card 2) 
(Card 1) 


'■ J J 0 U 0 5 fi U ? i v '■ •/ 1 1 P P f: i? \ 0 o 1 G G 0 0 D 0 0 0 0 D o 0 0 0 (ffi [f (3 B IS 0 B & 0 Q (] 0 0 i 0 ? 0 (? 0 0 U 0 {| 0 3 0 (J 0 CiO 00 0 0 0 0" 0 

1 1 ' 7 k ? II u nr r if ?. kv :*? h ^ Ji k ii * 37 » n <a * *: 4 . « 43 « <r At ai y a 5 ? m m m u k u a e i « » ;• vi?i .4 tj k jt :i tt to 

i n !,i n n ! 1 n 1 1 1 n 1 m i 1 n n m 1 1 1 m 11 1 m m 1 m 11 1 n m 1 1 1 1 1 1 m 1 1 m;: 1 1 m u 

: m ? ? V;. v ? 3 2 ? ? ? 2 z 7 23 a 2 2 2 5 2 2 3 2 2 * ? 2 ? 2 2 2 ? 2 1 2 2 2 2 2 2 2 2 2 2 2 2 2 1 2 2 2 2 2 2 2 1 2 2 2 2 2 2 2 2 2 2 2 ( 1-2 222222 

‘ • 4 K< 3 3 a M 3 3 3 3 3 S 3 3 3 3 3 4 3 3 * 3 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 * 2 3 3 3 3 3 3 3 3 33 3 3^3 3 3 3 3 3 i 3 

4 5 < 4 * < 4 < 5 -I 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4:4 4 4 4 4 4 4 4 [ 

; ■ 3 5 5*5,6 fi & S & 5 5 r . 5 ?■ ; J 6 5 5 5 5 6 5 5 5 5 5 5 5 5 0 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 S S Js 5 5 5 5 5 5 5 1 

i'll | 

I fi 7 r. fr :' >• & fi 6 6 fi 7 3 F * £ b £ e 7 F, f r ? M t C 6 r 6 6 ? 6 6 E 6 6 G 6 B 6 6 C 6 6 F 6 6 6 6 0 G 6 6 6 fi 6 G 6 6 6 6 6 S 6 6 G G 6'5 6 6 6 S 6 6 & 

8 7 7 7V7 / 7 > 7 7 V H 6 7 M U 7 7 7 7 7 7 7 7 Y? 7 M 7 7 i l 7 7 7 7 7 7 7 7 7 7 7 7 W 7 7 7 7 7 77 7 7 7 m 7 7 7 7 7 7 7 ! 7 7 7 7 7 7 7 7 | 

& 5 *■’? ■’ • - i* t 5 r : E r t t fl a 2 8 a f» U 8 B 8 li P S ft f 8 8 6 £$ 8 D 6 ft 6 6 fl T 8 & 6 8 8 8 5 0 0 3 5 i 8 8 8 8 8 8 8 3 3 8 B P-$ ? 8 8 C B 8 8 ! 


V 


■ 5V r 'j t 


: f . ■; C- 4 0 v F : ( g (| 5 9 ? S R R 9 (} 9 t) !5 9 n $ 3 0 S ? S 0 N g 9 5 g 9 9 9 55 a 5 9 3 9 » o 9 g 3 ? q 

; -w v > m :i L 2 :j r. h ■ ' ?: ^ 1 m fi (,,* , : - . ,? r ’& - ' ‘ j; , \ ,• 


This data set was obtained from a test on a graphite -epoxy tube. 
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The output listing for problem 2 follows: 

PROBLEM 2 COMPRESSION TEST 

COMPRESSION OR TENSION TEST ON TUBE SPECIMEN 
USING GERBER PLOTTER 

TUBE NUMBER 2 


FRAME 

5 

10 

IS 

20 

• 

75 

80 

85 

90 

95 

100 

105 

1 10 
115 

• 

150 

155 

160 

165 

170 

175 

THICKNESS = 


LUAD 

1 . 38b5E*J2 
5.23G8E+U2 
9.435dE*G2 
1 * 35 79E ♦03 

• • 

1.34636*03 

1.80596*03 

2.18l4t+03 

2.61546*03 

3.04336*03 

3.46246*03 

3.85376*03 

4.3034E+03 

4.74616*03 

• • 

3.06956*03 
3.51056*03 
3.93056*03 
4. 3392 E*03 
4. 74356*03 
5.0212E+03 
5.150 OE -02 


SG-1 

1.58416-05 
4. 1687E-05 
6. 83676-05 
9. 1712E-05 

• • 

9.6715E-05 

1.21736-04 

1. 45916-04 

1.6925E-04 

1.9176E-04 

2.192BE-04 

2. 41 79 E — 04 

2. 6847E-04 

2 . 9348E- 04 

• • 

1.9676E-04 
2. 2 1 7 8 E - 04 
2.46 79E-04 
2.70136-04 
2.9265E-04 
3.0932E-04 
OUTER 


SG-2 

3.83526-05 
1. L0056-04 
1.9U93E-04 
2.6597E-04 

« • • 

2. 6513 E— 04 

3 . 48 51 E— 04 

4.2021E-04 

4.9942E-04 
5. 8029 E- 04 

6*00336- 04 

7 • 3453 E -04 

8.2041 E-04 

9 ■ 0378E- 04 

• • • 

5- 8863 E-04 
6.72006-04 
7.5204E-04 
8. 3041 E-04 
9.05456-04 
9. 538 IE- 04 
DIAMETER * 


3* U 1506*00 
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STRESS-STRAIN POINTS FOR SG- 1 FRAMES 6* THROUGH 119 


FRAME 

STRESS 

STRAIN FRAME 

STRESS 

STRAIN frame 

STRESS 

STRAIN 1 

-RAH6 

STRESS 

STRAIN 


KSI 

MICRO 


KSI 

MICRO 

KSI 

MICRO 


KS! 

MI CRQ 



IN/IN 



IN/IN 


IN/IN 



IN/ IN 

64 

1.07u96+JO 

4*41096*01 

65 

1*09466*00 

4.66906*01 66 

1.29 67E+00 

5.065 9E +01 

67 

1. 4534E+00 

5.50276*01 

3d 

1 .62096+00 

5. 83626+01 

69 

1- 8031E+0Q 

6.67006+01 70 

1* 99986*00 

7.33 70E+01 

71 

2.18556+00 

7 .92066+01 

72 

2.35496+00 

8* 3375E + 01 

73 

2.5080E*00 

6.67108*01 74 

2.65546*00 

9.0879E+O1 

75 

2.8121E+00 

9*671 5E+01 

7 6 

3.00&9E+00 

1 *0338 E*02 

77 

3,22386*00 

L.0839E+02 78 

3. 4367E *00 

1.1339E+02 

79 

3*50986+00 

1 -20066+02 

80 

3. 76&4E +00 

1 * 21 73E+ 02 

81 

3.9231E+0O 

1.26738+02 82 

4.07 43E*00 

1.30066+02 

83 

4* 22906*00 

1.3590E+02 

84 

4*38386*00 

1-4257E+02 

65 

4.5456E+0Q 

1.45916+02 86 

4- 71 71E*O0 

1.5007E+02 

87 

4* 9U29E + OQ 

1*55086+02 

85 

5.0904E+UQ 

1* 584 1E + 02 

89 

5*2 744E*00 

1.64256+02 90 

5.45486*00 

1.6925E+02 

91 

5. 6350E+00 

l . 7342 E + 02 

92 

5 * 0 1JQ6 + U0 

1.75096*02 

93 

5.99206*00 

L.8342E+02 94 

6.16666*00 

1 . 88436+02 

95 

6. 34726+00 

1 .91 76 E +02 

9 6 

6 • 5201 E * UO 

1.9843E+02 

97 

6* 6951E + 00 

2.03436*02 98 

6.87536*00 

2.0644E+Q2 

99 

7 . Q484E+G0 

2.13446*02 

100 

7.22136*00 

2.1928E+02 

101 

7 *3 8526*00 

2 . 24286+02 102 

7.54926*00 

2-28 456+02 

103 

7. 7 1506+00 

2.3428E+02 

104 

7*877 Ot ♦ JO 

2*3 7628+ 02 

105 

8. 03?4E *00 

2. 41796+02 106 

8.21766*00 

2.45966+02 

107 

8.41246+00 

2 . 5346E+02 

108 

8*59e>4E + 00 

2 *5763 6 + 02 

109 

Q * 78 7 6E *00 

2-63466+02 110 

8.97536+00 

2.68476*02 

111 

9. 1 647E+00 

2.75976*02 

112 

9,34688+00 

2. 77646+02 

113 

9-53456*00 

2. 5181E *02 114 

9.7l66E*00 

2-88486+02 

115 

9. 89666+00 

2.9348E+02 

l 1 6 

1 « 00836 *01 

2.9681E+02 

117 

1 .02286*01 

3 *0098 E *0 2 118 

1 .03816+01 

3.043 26 +02 

119 

1- 0467 E +01 

3.0932E+02 

• 

• • 

• • 

• 

• • 

• • • 

• • 

• • 

• 

• • 

m m 


STRESS-STRAIN POINTS FOR SG- 9 FRAMES 64 THROUGH 119 


FRAME 

STRESS 

STRAIN 

FRAME 

STRESS 

STRAIN i 

FRAME 

STRESS 

strain 

FRAM6 

STRESS 

STRAIN 


KSI 

MICRO 


KSI 

MICRO 


KSI 

Ml CRO 


KSI 

MICRO 



IN/IN 



1 lx/ I N 



IN/IN 



IN/ I N 

64 

1 *07096 +00 

1.6675E+01 

65 

1 . 0 94 6fc *00 

1.63426*01 

66 

1 .29 o7 6+00 

2 .16776+01 

67 

1.4 534E + 00 

2.66806+01 

68 

1*62096+00 

3.33506+01 

69 

1.8031E+00 

3*91866+01 

7 U 

1.99986+00 

4.5022E+01 

71 

2. 1 855E+00 

5.00256*01 

72 

2.35<*9E*0U 

a.4194E+01 

73 

2.50806+00 

3.9l9ot*Lll 

74 

2.05346+00 

6.41 99E+01 

75 

2.8 1216+00 

6.83b7E*0l 

7 6 

3* 0069E + 0 J 

7.50376+01 

77 

3 .22386+00 

d . J J 4J 6 *0 1 

7b 

3.43676+00 

8.67106+01 

79 

3.6098E+O0 

9.25466*01 

80 

3. 76646 +00 

9.O715E+01 

61 

3.92316+00 

1 . Ul 72t *02 

«2 

4.07436+00 

l *07556+02 

83 

4.22906+00 

1 .11726*02 

34 

4. 38386+00 

1.1539E+02 

85 

4. 54566+00 

1.21 73E *02 

86 

4*717-16 + 00 

1 .25906 +02 

87 

4.9029E +00 

1 * 300 6 E *02 

88 

5.09046+00 

1.39246*02 

09 

5.27446+00 

1.4340E+02 

90 

5.4548E+00 

1 .4841E+02 

91 

5.6350 E+ 00 

1 .55916+02 

92 

5* 8 1006 *00 

1.60086+02 

93 

5.9920E+00 

1.66756+02 

94 

6.18686+00 

1.72596+02 

95 

6. 3 4 72£ + 00 

1*76756+02 

96 

6.52016+00 

1.82596*02 

97 

6*695 16+00 

1.8643E+02 

98 

6.8753E+00 

1.95936+02 

99 

7.04846*00 

2.00936*02 

• 

• • 

• • 

• 

♦ • 

• • 

• 

• • 

• • 

• 

• • 

• • 


STRESS-STRAIN POINTS FOR AVG OF GAGES 258-0 FRAMES 04 THROUGH 119 


FRAME 

STRESS 
KS I 

STRAIN 1 
MICRO 
IN/ IN 

FRAME 

STRESS 

KSI 

STRAIN 1 
MICRO 
IN/IN 

FRAME 

STRESS 

KSI 

STRAIN 
MICRO 
1 N/IN 

FRAME 

STRESS 

KSI 

STRAIN 

MICRO 

In/ in 

64 

1.07096+ OU 

9. 61596*01 

65 

1.0946E+00 

1.0ll6fc*J2 

66 

1 .29 6 7 £+00 

1 .1 6676*02 

67 

1.45346+00 

1.35346+02 

68 

1.&209E+00 

1.50636*02 

69 

l . 803 lt*00 

1.6925E+0 2 

70 

1.9998E+00 

1*07596+02 

71 

2. 1 055E+OO 

2 .03996+02 

72 

2.35496*00 

2.19006*02 

73 

2.5 OflOE+OO 

2*33736 +02 

74 

2.69546+00 

2 .47906+02 

75 

2.61216+00 

2.60966+02 

7 6 

3.00696 +00 

2. 8 J 706 + 02 

77 

3.2 236E+00 

2.99596+02 

78 

3. 43676 + 00 

3.19056*02 

79 

3.00986+00 

3 .34616*02 

80 

3. 76646 + 00 

3.49066*02 

81 

3. 92316+00 

3.6296E+J2 

82 

4.O743E+U0 

3. 77696+02 

83 

4.2290E+00 

3.9103E+02 

84 

4 .383 86 + JO 

4.06076*02 

85 

4.5496E+00 

4*2243 E ♦ 02 

06 

4.71716*00 

4.3627E+02 

87 

4.9 029 E + 00 

4. 55506+02 

88 

5.0904E+00 

4.73016*02 

89 

5.27446+00 

4.9J25E+J2 

90 

5.45486+00 

5. 06646+02 

91 

5. 63506+00 

5. 23076+02 

9 2 

5.81O0E+J0 

5.40276+02 

93 

5 *99 20E+00 

5.57506+02 

94 

6. ltiCibE + OQ 

5*73626+02 

95 

u* 34 726 + 00 

5 .90306+02 

96 

6.52016+00 

6*06976+02 

97 

6.69516+00 

6.2365£*u2 

98 

6.07536*00 

6 .40606+02 

99 

7.0484E + 00 

6.56446+02 

100 

7. 221 3E + JO 

6.73&7E+02 

101 

7.38526+00 

0*89 516 *U2 

102 

7*54926+00 

7.0480E+02 

103 

7. 71506*00 

7.2092 E+02 

104 

108 

7.87706*00 

7.36766 + 02 

105 

8.03 746*00 

7.51 7bE*02 

106 

b .21 76E + UQ 

7.6900E+02 

107 

6.41246+00 

7.88176*02 

8. 5964E+0G 

8.06246+02 

109 

0 * 78766 + OU 

0. 23466+02 

Llo 

8.97536*00 

6.40986*02 

111 

9. 164 7E + 00 

6.59596+02 

LI 2 

9.34&8E +00 
1.00836+01 

8.77946+02 

113 

9.53456*00 

0. 95176+02 

114 

9.7166E+00 

9.1266E+J2 

115 

9.09066*00 

9.29916+02 

116 

9.46866+02 

1 17 

1*02206+01 

9.5*376+02 

118 

1 -O381E + 01 

9.7578E+02 

119 

1. 04676+01 

9.81316+02 
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STRESS STRAIN TANGENT MODULUS 

1.0709E*03 9.6159E-05 l.O870E*0? 

X. 09*65*03 L « 01 16 E-04 1.08 67 6*07 

1.29 67E*Q3 1.18676-0* l-0858E*C7 

l.4534E*03 1 .3 534 E-04 1.08406*07 

1.62096*03 1.50636-04 1.08406*07 

1.80316*03 1.69256-04 1.0830E*07 

1.99986*03 1.87596-04 1.O82OE+07 

2. 18556*09 2.03996-04 1.08116*07 

2. 35496*03 2.19006-04 1.08036*07 




9. 71666*03 
9 • 890 6E* 03 
1.00036*04 
1.02286*04 
L • 03 8 IE * 04 
1.04676*04 


9. 12686-04 
9.299IE-04 
9.4686E-04 
9.59376-04 
9.7578E-04 
9.81316-04 


1.0427E*07 
1. 04 18E*07 
1.04096*07 
1. 0402E+07 
1.03936*07 
1 . 0390 E * 07 


STRESS-STRAIN POINTS FOR AVG QF GAGES 2 5 8 -a FRAMES 64 THROUGH 119 


FRAME STRESS STRAIN FRAME STRESS STRAIN J 

KS I MICRO KSI MICRO 

IN/IN IN/IN 

64 1.07096*00 9.61596*01 65 1.0946E+Q0 L.Uil6t*02 

68 1.62096*00 1.5063E*02 69 1.80316*00 l.o925E*Oi 

72 2.3 54 9E ♦ 00 2.19006*0? 73 2.50BQE*00 2.33736*02 

76 3.0069E+00 2.8070E*02 77 3.223BE+00 2.99596*02 

80 3. 7 6646*00 3-49066*02 81 3.9231E*00 3.62966+02 

S4 4.38386*00 4.06876*02 85 4.5496E+00 4.22436*02 

88 5.09046*00 4.7301E+0Z 69 5.2744E+00 4.90256+02 

92 5.01006*00 5.4027E*02 93 5.99206*00 5.p750E*Q2 

96 6. 5201E +00 6. U697E *0 2 97 6.69516*00 O.23&56+02 

100 7.22136*00 6.7367E*02 101 7.3852E+O0 6.89516*02 

104 7.8770E+00 7.3676E+02 105 8.0374E*00 7.5l7b£*U2 

106 8.6964E+00 8.0624E*02 109 8.7676E+00 0.234b6*u2 

112 9.34686*00 8.77946*02 t 1 3 9.5345E+00 6.9517E*02 

116 1 .Q0B3E *01 9.46866*02 117 1.0228E + 01 9.09376*02 


fc - 1. 04Z22e*07 FOR FRAMES 70 THROUGH 115 


RAMfc STRESS STRAIN FRAME STRESS STRAIN 

KSI MICRO KSI MICRO 

IN/IN IN/IN 

66 1.2967E + 00 1 • Id6?E*02 <>7 1.4524E+00 1.3534E*02 

70 1 .5998 t *00 L.8769E*02 71 2.18556*00 2.0399E+02 

74 2.65646+00 2.4790E+U2 75 2.81216+00 2.60966*02 

78 3 . 43676*00 3.19056+J2 79 3.6098E*00 3.34616*02 

bZ 4.0743E+00 3.7769E+U2 83 4.22906+00 3.9103E+O2 

86 4-71716*00 4.3827E+02 87 4.90296*00 4-5550E+02 

9U 5 .4540 E * OO 5.06646*02 91 5.63506+00 5.2387E*02 

94 6.16 68 E*00 5 .7362 E *02 95 6.34726*00 5.9O3OE*02 

yS o. 87536*00 6.40606*02 99 7.0484E*00 6.56446*02 

102 7.04926*00 7 .04 806*02 103 7.71506*00 7.2092E*02 

106 6 .2176E *00 7.6VOOE+OZ 107 d.4l24E*00 7*8817E*02 

110 tt_9753E*00 8.40986*02 111 9.16476*00 8.59596*02 

114 9.7 l6fat+00 9.1268E*02 115 9.89866*00 9.29916*02 

118 1 . 03 8 1 £ *01 9.75786*02 119 1.04676*01 9.61316*02 


NUMBER of plots - 3 

LENGTH OF A AXIS * 20 INCHES. 

LEN5TH OF V AXI S « 14 INCHES. 
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Problem 3 


This example illustrates the input, output, and plots obtained from the torsion test 
on two tubular specimens. The first tube was instrumented with three strain rosettes 
(each rosette contained three strain gages) and four LVDT’s. A stress -strain curve 
based on the average absolute values from all ±45° oriented strain gages and a torque - 
twist curve based on the LVDT outputs are plotted in figures 7 and 8, respectively. The 
second tube (number 4AL) was instrumented with two strain rosettes. A curve based on 
the sum of the absolute values from the ±45° gages (case 9) is plotted in figure 9. The 
data are expressed in U.S, Customary Units, A listing of the input data follows: 


/. 


A 


/ 


/■■"'ll 3 -“ *"• 

1 0. 


- 0041-5 

\ r \ 

! 

(Card 5) 

/ : * 1 iJ ! in 

(Card 4) 

y SUBE| JriUnSLK 4ftL 




1 <v 





1 > 

(Card 3) 
(Card 5) 

/ . l\Z\3 r y . 0 1 9 4 h- o fi !?, sr f! n 

7 ^ 

11. u !J f j 0 

1. 0 

. 00073 

! ” Tn 

/ lo \7[ 7 5 3 2 1 7 




4 N 

(Card 4) 
(Card 3) 

/ | |TUBt MUMK ■ \ 

\ > 



■* 

1 N 

/ 1 j 1 0 C r 1 1 20 14 1 

' 1 ! 

/PftfltJlJEM 3 T HRS* I OH TE5T 

0 


' 

a 

(Card 2) 


:1 . 
■/'sh 


fortpan statement 


s»JC O-jO. '0; I 0 0 U 0 

- ** 

0!D 0 !) u 0 si 0 0 fl 1 3 s 2 0 0 v'O 1 V'l t 0 0 3~0 0 0 !; 0 0 0 0 Q ii U ti C 0 P C S i; 

7 * * MM Ml IJ si .i :i 11 : ? n j.» n » M ji *1 r, y. a K „ r . , 


GOiJCfiJUuPGt 0 G f ft L i > 0 Co t 0 ,7 C Q D 5 

V u 'jV "'iniiimiimiiiimm , 1 Vi' , 1 , i m 7 m m n mm nnm n' mji’m mm 


3.3 3 3 3;*jJ J 3 « 3 3 4 3 1 3 3 3 3 3 4 3 3 4 3 3 3 3 3 S 3 3 3 3 3 3 3 3 3 3 7 3 3 3 3 3 3 3 3 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 J 3 3 3 3 s'l 3 3 3 J 3 3 3 

4j4 4 4 «j4 S J 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 1 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 * 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4'4 4 4 4 4 4 4 4 

SjS 5 S 3|Sj( 5 5 3 5 5 5 3 5 5 5 5 t 5 5 6 5 5 5 5 5 5 5 5 5 5 5 5 5 5 S 5 5 5 5 5 5 5 5 5 5 5 5 5 5 S 3 5 a 5 5 5 1 5 5 5 5 5 S a 5 5 5 b 5 f/j 3 S a i a 5 5 

*j* 6 7 5 ;f 5 s ‘ * 5 M 6 i S 7 6 8 6 6 6 C 5 5 B 5 3 C 4 S S 6 C 3 6 C 6 t D 6 8 c E 5 6 S 5 0 0 8 $ 6 6 S 5 4 6 S 8 6 4 G S 6 6 8 5 5 6 0^6 i S 6 8 S 6 I 

7 j S 1 1 7 j J j 7 J ' 7 7 7 7 J I 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 J 7 7 7 7 7 7 7 7 7 7 7 7 7 7 ?;// 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 / 7 ?'j 7 7 J 7 7 | ; 
i S-3 « «) *33 33«.3l3 9 3S-»aB#!msi«i3*e»Si68iStS?:M4 3?»$J43J3 3 3 3 0 JS3 < 1si*i*3Jt , *flt3*»i 

5 3 3 3 i 3 3 f 3 1 i 3 f 9 3 a 9 i 3 S 3 5 3 « 5 j 3 i. 4 ■!«. • ■ - : . : ■. a s i - - -- 0 1 . - • S S 3 s ' 9 9 9 i : S L - 


i N (Card 1) 


l I 

i ■ j :• s S'j a i <! -j i ; 4 • 

i i . ; i s • ■ ; i ■) : 

. PpQr.W 


^5?B339*8S3S2S3j3b5 0 • :■ 

■ . -c j . jj a j<i ' 5 : u x 0 ! -i 1 : -i •.* <j y , 


The first data set was obtained from a test on a graphite -epoxy tube whereas the second 
data set was obtained from a test on an aluminum tube. 
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The output listing for problem 3 follows: 


RRJBtEM 3 TORSIUN test 

TORSION TcoT UN TU0 l- SPECIMEN 
USING bcRUeR PLUTTtK 

rUb£ NUMliEk 4 


5 

10 

15 

2J 

c5 

j3 


LuhO 

b .2414 l+ 0 l 
i - 41286+02 
3* 730<tt: *02 
5 • 95 4d t *0 2 
8. 35 94 c * 0 2 
1 .08 48c +0 3 


5G-1 

l .4 l 74E-05 
4 . JU2Cc-05 
/.a 7e5L-05 


SC-2 

-3 . 3 J 50c~06 
-8. 3 375E-06 
-1 . 75U9 l-05 


1.2089L-04 - 3 .00 1 5c- 05 
l. 6925 — 04 -4.U054F-O5 
2.2094E-O4 -5.33602-05 


• ♦ 




• • 


160 
1<>5 
170 
i 76 
1 iO 
id 5 

THlCKNcSS = 


2.1452E+03 4.3,772c-U4 -1.0589E-04 

2 *4<« 7t» u. ♦ 05 5.bl42E-J4 -1.158SE-04 


£. 75571+03 
3. J/23t+05 
j. Jdfl'lc^Oi 
5 .0267 l +0 i 

5 .3 700L-J2 


S.olLlc-U* -l*ii4JE-04 
6 . 2 7 d l E- J4 -1.459iii-04 
6. 42856-04 -L.5B4L6-04 
7.4JJ7E-J4 -I.7CCEE-U4 
JUTER DIAMETER = 


3.0 14 4 E +00 


AVG STRlo^-STk AIN RJlNTS P Jk FRAMES 7 THROUGH 75 


FRAME 

jTRlSS 

$ 1 r\^ L N FRAME 

STRESS 


KS 1 

MICRO 

KSI 



I N/ IN 



STRAIN 

M1CRL 

tN/IN 


7 

1 . >4Ult- 01 

4. 3074c +31 

a 

2.92 42c>ul 

7.64436+0 1 

15 

5. JO 04t- 0 1 

i .uJ i4t.+ J2 

19 

7.4490E-di 

2. J734l*02 

2 3 

■jf, 79&7C.-0 l 

3. Lt>2 U+J2 

27 

1 .2 5 75c *00 

•t • 02 9 3l + 02 

Jl 

I . 64 2 JE *00 

4 • 9 2 *♦ 7 c + J2 

35 

l.d J42ii + 0<J 

V * 8640= + d2 

39 

2.13 7 ?c + 0 J 

b . 82 2 tit ♦ 02 

4 J 

2 « 4 *i 1 ?L ♦ 00 

7 . ttO 44 l + 02 

47 

2 . 7 5 3 Ot * 00 

rt. 7 8 4 4 c + J 2 

51 

i.Ot> l>5l *• 00 

9 • 76 82 E + Oil 

55 

3 . > 7 4 a c. + u J 

i . j78jc+ J i 

54 

J. /J2dc+00 

i . U 1 1L *03 

o3 

JO 7 3c + l)0 

1 .28542+03 

6 t 


i.,*9 2K + UJ 

71 

4 . a 4 4 8 L + J 0 

1 .4 4 6b c + J J 

75 

4 , d-t 4 7E *-00 

l . b3 8Ji_ + J3 


8 

l* 7707E-01 

5. 752 BE +0 I 

i* 

3 . 4t24 E- U l 

1 .1200E+02 

16 

5, 57 UE-UL 

1. 8065E+02 

2 J 

7. 9494 t-0 I 

2.5741E+02 

24 

1. 06J2E+U0 

3. 3 96 l L + 0 2 

23 

1.326?£+00 

4.2549L +02 

32 

1.61 L3E+0U 

5 . 1 6 J 7E +0 2 

Jo 

1. 91 12E+00 

o . 1030C+02 

4d 

2. 21 C4c+G0 

7.056 31 +0 2 

44 

2.5165 fc+CC 

8.0374E +02 

48 

2. £3 17E + CC 

9. 0406E+02 

52 

3. 145 JE* CO 

i . 0JJ3E+U3 

56 

J. 4 64 5 E + 0 Q 

1 . 1 036£ +03 

to 

3. /8 74 t + 00 

1.2U67E +03 

64 

4. 1141E+0J 

1. 3L26E+03 

6 d 

4. 4440C+CO 

1 .4190E+03 

72 

4. 75 63 r + JO 

U5L J8T+03 


O = 3.14202L+U6 PuW FRAMtS 7 THROUGH o5 


A Me 

STRESS 

KSt 

STRAIN FRAME 
MICRO 
I N/IN 

STRESS 

RSI 

4 

2. U68E-01 

6 * 9 7 5 6 E +01 

LO 

2 . 5 7B2 1-0 1 

l 3 

3.9624E-U 1 

1.2895E+02 

14 

4.5006E-01 

1 7 

6, 1923E-01 

1.9927E+02 

LS 

6.8J75E-0L 

21 

8.6530E-01 

2. 7930E+02 

22 

9* 3451E-01 

25 

1 . 1267E+Q0 

3.6101E+U2 

26 

L . 1 998 E + 00 

29 

U3998E + 00 

4.4828E+02 

30 

1.4b8 9E + 00 

33 

1.6881E+00 

5 . 3888E +02 

34 

L.7688E+00 

37 

i *9880 £+00 

6.3J65E+02 

38 

2.0571E+00 

41 

2.29L7E+00 

7 * 3064 1 +02 

42 

2 • 3 u8 6 E + 00 

46 

2 « 5 4 93 T. +00 

8 . 2 9 02l +0 2 

46 

2.6 72 3E + 00 

49 

2.9106l* + 00 

4. 2880c +02 

50 

2.9915E+00 

53 

J.2l83t+00 

1 .0280c +0 3 

54 

J. 3068E + 00 

57 

3.5413E+00 

L. 1289E+03 

58 

3.622 OE +00 

61 

3*8719 £ + 00 

1.2337c +03 

62 

3 . 960 3E +00 

65 

4. 1986E+00 

1 .3387E+03 

66 

4.2872E+00 

69 

4.5J7LE+00 

1.4454E+03 

70 

4.62186+00 

73 

4 « 806 3 E +00 

1.5266E+03 

74 

4. 8332E+00 


STRAIN 
MICRO 
IN/ IN 

8. 3Q97E+01 
1.4507E+02 

2. 18996*02 
2.9653E+02 

3. 8U74E+02 
4.6995E+02 
5.6334E+02 
6.5B38e+a2 
7.5537E+02 
8.5348E+02 
9 . 538 1 E + 02 
1.05L9E+03 
1. L547E+03 
1.2590E+O3 
1. 3654E+03 
U4713E+03 
1- 5350E+03 
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TORQUt TWIST POINTS. FRAMES 


1 THROUGH 75 


FRAME TOftOUc TnIST FKaMfc TLRCUt TWIST FRAME TORQUE TuIST FRAMc' TORQUE TWIST 



IN-L8S 

0 E o/ IN 


IN-L0S 

UtO / J N 

7 

1. 1426E+02 

1 .d359c-0J 

3 

1.3137E+02 

2 . 2 1 3 3E -0 3 

LI 

2 . 16956+02 

.>. 72226-03 

12 

2* 56 38E + 02 

4.23606-03 

15 

3. 73d4E+02 

6.187 Ot- 03 

l 6 

4*13776+02 

6. d387fc-03 

11 

5*56406+02 

9.0248<_-03 

20 

5* S34BE + 02 

9. 7439E-03 

23 

7.41HIE+02 

1.20J4E-J2 

24 

7. SO JOE +02 

1 .2d33E-02 

27 

9.32 93E *02 

1.5273^-02 

28 

9. 8427E+G2 

1 . 6O96C-0 2 

31 

l. 1440E+03 

1 • 86 62l- 02 

32 

1. L954E+0J 

1.95 /86-02 

35 

1.3j 086 + 03 

2.2 2306-02 

36 

1.4179E+G3 

2. 313 IE -02 

39 

l.3d6lctJ3 

2*37ddfc-0Z 

40 

l • 640 3E +03 

2.67161-02 

43 

l.dl 15E+03 

2 . 95U3E-02 

44 

1 . 86d5 E + 03 

3.04J1F-G2 

47 

2.U425E+03 

3.32 7 2c- 02 

4d 

2. 1053E + C3 

3.4215f-02 

51 

2.2 7656 + 03 

3. 7095E-O2 

52 

2. 3335E+C3 

3.805 Lf-OZ 

55 

2 . 5 J 756+03 

4. 0 8 94 u - 02 

56 

2.5703E+03 

4. 1865E-02 

59 

2 • 74 7 it + 03 

'♦ *4844t-02 

6 J 

2*8059E+03 

4.58006-02 

63 

2*99336+03 

4 . ddo^c- 02 

64 

3. U5Z3E+OJ 

4. 9906E-J2 

67 

J.2377E+03 

5.2999E-J2 

6iJ 

3.2976E+C3 

5 ♦ 4Gri2I -02 

71 

3.48 31 E + 0'3 

5. 70 70l- J2 

72 

3. 5267E+C3 

5 . 77 98E-0 2 

75 

3.3943E+J3 

3.d776E-J2 






IN-LdS 

[)EG/1N 


IN-LBS 

0 EG/ IN 

9 

1.5705E+02 

2. 7207E-03 

10 

1 • 9 128 E + 02 

3. 1 796E-03 

13 

2.9397E+02 

4.9076E-03 

14 

3.3 390 E * 02 

5* 500 9E-03 

1 7 

4.59412*02 

7*5 224E -0 3 

13 

5 . 0505 c + 02 

8.282 JE-03 

21 

6*41 97 F* 02 

1. 05046-02 

22 

6 ■ 933 2 6 + 02 

1. L232E-02 

^5 

d.3594E+02 

L .36386-02 

2b 

8.9014E+02 

1.4434E-02 

29 

1.0335E+03 

1.6933C-02 

30 

1 • 0099 £+0 3 

1.77936-02 

33 

1.25246+03 

2.0439E-02 

34 

1.3123E+03 

2. 1313E-02 

37 

l .47496+03 

2.4C32E-02 

33 

1.5262E+03 

2.49L6E-02 

41 

1.70026+03 

2* 7644E-02 

42 

1. 7573F+03 

2.85716-OZ 

45 

1 . 92 64 1 + 0 3 

3* 1 4 3 It -02 

46 

1.98266+03 

3.2293E-02 

49 

2. L5V5E+03 

3.5209E-02 

50 

2.2L94E+03 

3.6123E-02 

53 

^ . 38 77 F+0 3 

3.6972E-02 

54 

2.4533E+03 

3*99595-02 

37 

2.6273E+03 

4.2874E-U2 

53 

2.6872E+03 

4*38 ?2E-02 

6 1 

2. 8726 iT + 03 

■t.6364E-02 

u2 

2 . 9382t*03 

4. /372E-02 

o 3 

3. L L 31 6 + 03 

5 ♦ 090 16-02 

6 6 

3. 1807F+03 

5. 19506-02 

69 

3*366 Lc + 03 

5.5 11 76-02 

70 

3.4289E+J3 

D. 61417-02 

73 

3 ■ 5o 381 + 0 J 

5. 83556-02 

74 

3 .58586+0 i 

5.870 85-02 


SLOPE J F CURVE = 6.L3d?dc+13 FOR FRAMES 7 THROUGH 65 


TURF *'UMHr R 4AI 


FR A*>F 

L'lAO 

SO-1 

s t- ? 

10 

3. 531 IF + 0? 

- 004 : r 

3. 

?0 

1. 173*r 5 

- 1 . < 0 3 5 r - )4 

-1 .3335 F-3 5 

3 ) 

1.961 ? t +n» 

769F-04 

0 

1 

40 

2.9 3 r \F +0 3 

- ‘3 . 59 4 5 C - [)4 

-■*.5 M 2P-C5 

50 

-.8 5 76 F’+ 33 

-7. -2 ? 7 c -*}4 

-3 . J 31 -05 

60 

4 . RS30T +03 

-3.07 *04 

-A. 6 3] 7F -05 

7J 

5. 339 7! +0 3 

- 1 . U5 0? 

-4 .835 7F-05 


• 

• 9 • • 

• » 

220 

1 . 692o c+04 -3.3225L-03 

1. 7842E-04 

z>0 

1.6735E+J4 -2.767^6-03 

6 .35326 -04 

2+U 

1.2ci7oE*04 -L.lddlE-JjJ 

3.1 2666-04 

250 

i • 1 2 40c+04 -9.5631E-04 

3 -66 8 5E-05 


THICKNcSS " 4.90 JOl-02 UlITEP. ItA^ETER = 3.000 Gt+OO 
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AYG bTRcSS-ST*A H\| PulNTS FJR FRAMES 4 THROUGH 220 


FRAME 

aTkdSS 
KS J 

STRAIN hRAMc 
rl 1 CKO 
l N/JN 

STRESS 

KSI 

STRAIN FRAME 

micro 

IN/ LN 

4 

5.ji92t -02 

<♦ . 5 it *0 l 

5 

1,41 iOE-Ol 

5, 7520E+O1 

6 

tt 

2 . 2t> 82c- 0 l 

1.15476*02 

9 

3.5539E-CI 

1.39246*02 

LO 

12 

6.65 39E-JI 

2.22 1 9c *02 

Li 

8.0966E-OI 

2.5 179l*02 

L 4 

lu 

l . it* 97c* JO 

3 . 4i 1 2c * 02 

17 

1 . 34 1 L c *00 

3-7560E+02 

18 

2U 

1 . U0 3 46 + Oo 

<♦.68 ir>t*02 

21 

1. 76 S 7 E *00 

4 . 96 9 IE +0 2 

22 

24 

2«2d4 OE + 00 

5 . d9bb2+02 

25 

2. 22 6b£*CG 

6.215oc*02 

26 

23 

2. 71 25l + JO 

7. 1619c *02 

29 

2 * 7982 E *C0 

7.48706*02 

3 0 

32 

3 . 53E + 00 

6, 504 Ji_ *02 

3 j 

3. 4839L*Q0 

6. 05066*02 

34 

36 

3 * dO 96t *uO 

■i.o7 V9ui*0 2 

37 

j, 26 46 t ♦ C J 

1.02226*03 

38 

40 

4 . Jd 3 ‘it! t(JO 

l.L23t>c*di 

4 L 

4- 56 53 E+CC 

1 . 1002003 

42 

<+4 

4.90396+00 

l . 26362*03 

4 5 

5. 026 7£*00 

1 . 29906*0 J 

46 

4 3 

3.4-* 92 L * OJ 

l .40 651*0 3 

49 

3. 7 1 24E *C0 

1 .4 J90L *0 3 

50 

52 

6 • 1 *t 1 Oc + 00 

1 . 54 70c + 0 3 

63 

6. 226 7£ *00 

1.5825F *03 

64 

56 

6. o0 53c * 00 

l.od92c*iJi 

6/ 

6.78 391: *00 

l . 72466*03 

68 

60 

7.2552l+00 

1 . 6 3 0 L c * 0 3 

Ci 

/. 34 10E *00 

1. 86 726 *0 T 

62 

64 

7./64-3L *00 

L.47 i‘Jc*0 3 

c5 

7. 85 5 J E *00 

2.01 J2F *03 

66 

6 a 

3 . 4 6 3 3 1 + 0 J 

2.12 Otic* 0 J 

09 

8.4553E+C0 

2. 1556c *u 3 

70 

72 

o. ->2 oo p . + J*j 

2 • 2l> 4 Jd *0 3 

73 

9.09806*00 

2 . 2994r *0 i 

74 

76 

9 » ** a 3 H c + 0 U 

2 .40 7 1 1 *0 i 

77 

9. 7 8 3 7 L * CO 

2.4458L + 63 

7a 

dj 

l .Jt 7 uL * 0 1 

2.66 4oii * J J 

d 1 

1. 0384E*Cl 

2.592 IE *0 3 

82 

34 

1 . *>0 36c *01 

2-^993l+J J 

06 

1. 094 IE *0 1 

2. 7J51E+03 

86 


STRESS 

KS] 

l.u LOE'Ol 
5 • 26 8 L£ - 0 1 
9.9023E-O1 

1 .4268E+00 
1 . 8982 d* 00 

2 ,54 1 GE *0 0 
2.92o7[-+00 
3. 4411000 
4 . 08 30E* 00 
4.7267000 
8 . 1553E+ 00 
5.7552c*00 
6.3981000 
6.9124E*00 
7 . 5 1 24E *-00 
8. 1123E+G0 
a . 7 L 24F *0U 
9. U24li-00 
9. 86 45l* 00 
1 ,047Ct*0l 
1.10276*01 


STRAIN FR\*E 

MLCRU 

IN/ I N 

STRESS 

KSI 

7 .4 20iE *0 1 

7 

L.4L10E-01 

1 • 663 JC *0 2 

11 

5.2601E-O1 

2.8L33E*02 

15 

9.9023E-O1 

4.08952*02 

19 

1.68396*00 

5.2 568L*0 2 

23 

2. 13546*00 

6.5157c*02 

27 

2*49d 2E*00 

7.84976*02 

31 

3. 10396*00 

9.2 004E *0 2 

35 

3.6 12 5E*00 

l . 0 56 8E *0 3 

39 

4. 2 124E *00 

1. L 952t *0 3 

43 

4. 7695 E*00 

I . 3 3 3 2t + 0 3 

47 

5. 241 OE *00 

1.47622*03 

51 

5. 041 OE *0U 

L .6 10 3C *03 

66 

6 . 4039 E *00 

1 . 76C0E *03 

59 

6.9982E*00 

l .90ldt*03 

63 

7.603 8 E*00 

2. 046 OH *03 

67 

0. L93LE+00 

2.19312+03 

71 

0. 798 2 E*00 

2 . 3 35 3 1: *0 3 

75 

9.390 1E*00 

2.482 9“*L)3 

79 

U 0041c *01 

2 . 62 75E *0 3 

a 3 

1 .06B4E+01 

2.77266*03 

07 

1 . 1284E + 0 1 


STRAIN 
MICRO 
IN/ IN 

9. 3797E+0L 
1 , 930 IE*02 
3. 12656*02 
4. 3055E* 02 
5 . 5 73 6E + 02 
6. 845 1 E*02 
8. 1 032E*Q2 
9.5339E+02 
1.09146+03 
U23L4E + Q3 
I. 3694E+U3 
1.51166*03 
L. 6 54 2E* 03 
1. 796 7 E +03 
l. 9331E+J3 
2.083 LC + 03 
2. 22 936*03 
2. 37L6E+03 
2.5192E+J3 
2. 662 IE * 03 
2. 91016*03 


• • • • • • 


10b 2.i50 3t+GL 

192 2.3384 c *01 

196 2. 4 J 5*ti: *01 

2JJ 2.422ut>Ji 
2JH 2.*Ou9c + 0l 
20ci 2.4ob^c*0l 
212 2.4u2/c*Ul 

2 1 1> ^.'*9 9 7d*Jl 

220 2.6-i66[L + 01 


0 . 64 1j 8E * 03 1 89 
ti.al'jJc+0 3 193 
6.9d22=*03 L 97 

1 m 1 5 Joe * 03 20 L 
/. j2 1 U*U3 206 
7.4d0i£*U3 2v,9 
7 .00 C4d ♦ J 3 2 1 j 
/. twiVE+Oi 2 17 
£3 . 0 3 36c *03 


• • 

• • 

• 

• » 

• + 

• 

• • 

2. 3503E *C1 

6 . 68 791. +0 3 

190 

2. 36266+01 

6.73086*03 

191 

2.38416+01 

2. 39696*0 1 

6. 05636*0 3 

194 

2.39266*01 

6.09896*03 

195 

2 .396 9E + 01 

2.414LE*0I 

7.O266C+03 

190 

2.4ia4h*Ul 

7 . 06476+0 3 

199 

2. 435 6E +0 1 

2.43S7E+C1 

7. 1 9 326*0 3 

202 

2.4569L*U1 

7. 2 35 3E +0 3 

203 

2 .439 7 E +0 1 

2. 46oSt+Cl 

7 • 36 L 6£ *03 

206 

2.4655E+0L 

7.40i3G*03 

207 

2. 4 74 1 E+0 L 

2.48696 *01 

7.53296+03 

2 10 

2.4027E*O1 

7. 5 730E+O3 

211 

2.4069E+O L 

2. 48696+01 

7, 70001 *0 3 

214 

2-49976*01 

7.75L8E+03 

215 

2.49976*01 

2. 5109E+01 

7 . 89 i 5E *0 i 

2L3 

2.50B4t*Ol 

7.93051+03 

219 

2.52122*01 






6. 77046*03 
6. 941 3E*03 
7. 1090E*03 
7.27566*03 
7.44616*03 
7.61926*03 
7. 7968E*Q3 
7. 9894 £ *03 


NJMBER of PLulS 
LENGTH Jr X AXIS 
LENGTH uF 7 AX I S 


3 

*0 INCHES 
14 IpNCHcS 


Langley Research Center, 

National Aeronautics and Space Administration, 
Hampton, Va., September 18, 1974. 
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APPENDIX A 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 

The International System of Units (SI) was adopted by the Eleventh General 
Conference on Weights and Measures in 1960 (ref. 1). Conversion factors for the units 
used herein are given in the following table: 


Physical 

quantity 

U.S. Customary 
Unit 

Conversion 

factor 

(a) 

SI Unit 
(b) 

Length 

in. 

0.0254 

meters (m) 

Load 

lbf 

4.448 

newtons (N) 

Modulus , 
stress 

psi = lbf/in 2 

6895 

newtons per meter 2 (N/m 2 ) 


a Multiply value given in U.S. Customary Units by conversion factor to obtain 
equivalent value in SI Units. 

^Prefixes to indicate multiples of units are as follows: 


Prefix 

Multiple 

micro (p) 

10 -6 

kilo (k) 

1— i 
O 
CO 

mega (M) 

10 6 

giga (G) 

10 9 

tera (T) 

1012 
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APPENDIX B 


SOURCE PROGRAM LISTING 


A listing of the Source Program BECKPLT is presented in this appendix. A 
description of the significant FORTRAN variables is given in appendix C. The matrix 
inversion subroutine MATINV used in this program is described in appendix D. 

PROGRAM rtECKPLTl 1NPUT = 20L, OUTPUT* 20 l, TAPE 5= INPUT, TAPE30*201) 

PROGRAM TO ASSIST IN THE REDUCTION OF BECKMAN DATA 

DIMENSION TITLEI(4),T1TLE2(4)»AC2(99) .IGAGESI4) »DTA(415»20)» 

IX14U0), Y ( 400 1 » TEMPX14), TENPYI4) TllkltTC 

COMMON GMGVfc, I PLOTS. NUPLT t 18 ) , XLN .YMULT.XMULT .1UNITS 

C 

c READ ♦ PRINT TITLE. 

READ ( 3,11 (TITLEII J),J=1»4) 

1 FORMAT 1 4A10) 

I F ( EUF , 5 ) 1000 , 2 

2 PRINT 3, < TI ILcU J ) » J =1,4) 

3 FORMaT UH1 4AIU // ) 

L REaD PLOTTER, PEN, INPUT TYPE, SPECIMEN TYPE, TEST TYPE, PLOT INFO, UNITS 

READIS,4) IPLT.IPcN.INPT.ISPEC^TtST.IXSlZE.lYSlZE.IOMOVE.IUNITS 

4 FGRMAT(9t4) 

IF( IYSIZE.EQ.O) IYSIZF=9 
IF! 10MOVE.EQ.0 ) IOMOVE=2 
I X= l XSI Z E 

XLN = I X SIZE *■ 0.0 
YLN = I YS l ZE *■ 0.0 
OMGYc.= luMOVE+O.O 

IFUTtST.EQ.O .AND. ISPEC. EQ.O) PRINT 41 
A1 FORMa I l 1H .^COMPRESSION OR TENSION TEST ON ELAT SPECIMEN*) 

IE( ITLST.EU.O .AND. ISPEC.EQ.l) PRINT 42 

FORMA T 1 1 H , *COMPRE SS I CN DR TENSION TEST ON TUBE SPECIMEN*) 

IFl I TEST .EQ.l .AND. ISPEC.EQ.l) PRINT 93 

43 FORMAT! 1H ,*T0RS10N TEST UN TUBE SPECIMEN*) 

IFIITlST.NE.O .AND. ITEST.NE.l) PRINT 44.ITEST 

44 FuRMAT ( IH ,*INVALID TEST NUMBER, ITcST = *,I4) 

IFUTiiST .Nc.O .AND. ITEST.NE.l) STOP 

1FI I SPEC .Nc.O .ANO. ISPEC.NE.il PRINT 45, ISPEC 
i,5 FuRMA TI 1 H ,*INVALID SPECIMEN NUMBER, ISPEC = *,I4) 

I F ( I SPcC . NE * D .AND. ISPEC. Nt.l) STOP 
IF! I TEST. EG. 1 -AND. ISPEC. EQ.O) PRINT 451 
451 F0RMAK1H ,*LANNOT HANDLE TCRSION TEST ON FLAT SPECIMEN*) 

1FI ITcSI • EU • 1 .AND. I SPEC. EQ.O) STOP 
I F { IP LT . EQ.O .AND. IPEN.EQ.O) PRINT 46 
^6 FOKMmTUH , *US ING CALCCMP PLOTTER WITH BALLPOINT PEN*///) 

1FI IP LT. cQ.O .AND. IPEN.EQ.il PRINT 47 
47 FORMA 1 1 1H f * US 1 NG CALCOMP PLOTTER WITH LEROY PEN*///) 

if; iplt.ne.oi print 4a 
4B FORMAT! 1H ,*USING GERBER PLOTTER*///) 

0 

c IN1 T IaLIZc PLOTTING ROUTINE. 

IF(IPlT ,EJ. 0) CALL CALCOMP 
IF( IPlT.Nc.O) CALL GERBER 

1 F ( I P LT . cU. 0 .AND. IPEN .NE. 0) CALL LEROY 
KOUNT =0 
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C 

C PRUT TITLE ON GRAPH. 

Call calpltk. , 1 . 2 , - 3 ) 

CALL NUT ATe I -3 . ,0.,.25,TITLE1(1),90.,40) 

c 

C READ «■ PRINT TITLE Cf- RUN. 

55 READ 15,1) ITITLE2IJJ), JJ =1,4) 

I Hit UP, 5 ) 1000, 1100 
1100 PRINT 37, {TITLED JJ) ,JJ=1,4) 

37 FuRMaTI IH , 4A l 0 ///) 

ChlL NUTATtlG. ,-i. 1, .20, TI TLE2I 1) ,0.,40) 

C 

C RcAD RUN INFO. 

1F(1Tc3T .tQ. 0) R EA D ( 5 , 6 ) J AC , NRL , NR U, IN Da VG , 1 GA GES , ROD , NRl SLP , 
INRUSlP,NGPLT,MORE 
o FDRMAT(4 13,4 12,313,1812, 14X, I 1) 

IFIlTEST.EQ.l) READ! 5 , 61 ) J AC, NRl. , NRU.NROS, I N DA VC , MOD, NRL SLP , 

1 NRU SI P »N OPLT , MORE 
61 FURMATl 3 13, 1812, 19X, I 1 ) 

UP = NKU-NRL+1 
1FIJAC.GT.20) PRINT 75,JAC 

75 FUR MAT ( 1H0 »* TOO MANY CHANNELS, JAC = *,14) 

I F ( JmC.GT.20) GU TO 1000 

c 

C ktAU IN DATA. 

IFUNPT .CO. 0) CALL ROT A PE ( JAC , NRL , NRU , MO , CT A ) 

I F ( IN PT .EL). 1) CALL RDC ARC ( JAC, NRL ,NRU »MD , CTA ) 

IFUNPT . EQ. 2) CALL R DC E L L ( J AC , NRL , N RU , M D , OT A ) 

KNTRUN=0 

C 

C READ IN SECTluNAL DATA. 

1FUTEST .EQ. 0) RtAD<5,8) THICK, WIDTH, SC FAC, RMAXSN 
a FORMAT l 4F 12.4) 

I F ( I TEST .EQ. 1) READ (5,81) THICK , W I DTH ,R LENGTH, AL ENGTH, SCFAC , 
1RMAXSN 

81 FURMaT (6F12.4) 

IF(SCFAC.EQ.O) SCFAC=1. 

82 IF(NP.LT.l) PRINT 76, NRL, NRU 

76 FORMA I ( 1H0 ,* LAST FRAME TC BE PLOTTED LESS THAN FIRST FRAME*, /, 

1* FIRST FRAME = *, 13,6X, *LAST FRAME = *,13, / , * WILL GO TO NEXT 
2PLUT* ) 

lF(NP.LT.l) GO TO 50 
C 

G USE TYPE UF SPECIMEN *■ SECTIONAL DATA TO COMPUTE STRESSES 
C AND STURc IN Y ARRAY. 

CALL LJAuYITHICK, WIDTH, SCFAC, ISPEC,ITEST, NRL, NRU, Y,DT A, FAC, YMULT, 

I IUNI T S ,K NT RUN) 

IFIMOO.NE.l) GO TO 751 

IF{ MOD.EQ.l .AND. NRLSLP .EQ. 0) NRLSLP=NRL 
IFIMUD.tQ.l .AND. NRUSLP.EC.O) NRUSLP=NRU 
IFINKLSLP.GT.NRUSLP) PRINT 613, NRLSLP , N RUSL P 
613 FORMATUHO,* FIRST FRAME FCR SLOPE GREATER THAN LAST FRAME FOR SLG 
lPfc.*,/,* FIRST FRAME = *,I4,6X,* LAST FRAME = *,I4»/»* WILL NOT FI 
2ND SLOPE.*) 

I F(NRlSlP.GT.NRUSLP) mqd=d 
IF ( MOD. t Q.U ) GO TO 751 

1FINRLSLP.LT.NRL .OR. NR LS LP .GT.NRU) PRINT 611, NRLSLP 
6ll FORMhTIIHO,* FIRST FRAME FCR SLOPE INVALID.*, / , * FIRST FRAME = 
l*, 14, / , * WILL USE FIRST FRAME PLOTTED.*) 

I FI NRLSL P.LT.NRL .OR. NRLSLP .GT.NRU ) NRLSLP = NRL 
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1FINkUSLP.LT. NRL -OR. NRUSLP .GT.NRU! PRINT 6I2.NRUSLP 
612 FORMA T l IHO »* LAST FRAME FOP SLOPE INVALID.*, / , * LAST FRAME 

114,/* vJILL USE LAST FRAME PLCTTED.*I 
I F l NI'.OSL P ♦ LT . NRL -OR. NRUSLP .GT.NRU J NRUSIP-NRU 
751 CONTINUE 


= *, 


SET A * V SCALES. 

XMULT = 1.0 

IFIRMAXSN .LT. .1) XMULT=IC.**6 
TEMPXU) =0. 

TEMPX (2) =RMAXSN* XMULT 
IFlISPtC.EU.l .«ND. ITEST. EG. I .AND. 
1 ( 2 ! 

TcMPY(l)=D. 

T cMPY ( 2 ) =Y l NP ) 

DO 7t> 1 1 = 1 , NP 

761 IF( Y( II.GT.TEMPYI2) J T EMPY I 2 ) = Y( I 1 


INDAVG.Nt.O) TEMPX12 I =2*TEMPX 


1 PLOT S-= l 

I F ( I T EST .EQ. 0 -ANU. INOAVG .EO. 01 IPL0TS=JAC-1 
I F( (TEST .tQ. 0 .AND. INOAVG .£U. 1) IPLOTS=JAC 

iFMSPCC .EO. I -AND. ITEST .EG. I .AND. INOAVG .EQ. 0) IPLOTS = 

1 3*NR0 S 

IFIISPcC .EG. 1 -AND. ITEST .EQ. I .ANO., INOAVG .EQ. 3) IPLOTS = 
13*NR0S*1 

I F ( IX.EQ.U .AND. I TEST .EQ.O I XLN=OMOV E* JAC 

IFIlX.LQ.O .ANO. ITEST. EQ.l) XLN=uM0VE*(3.*NRUS*l .) +NRQS 

IXSIZc-XLN 

XXLN = X LN-OMUVt* I IP LUTS-1 ) 

1FIXXLN .LT. UMUVt) XXLN=CMCVE 

IF( ITEST. EQ.l . mNO. ( INOAVG.EQ.O.GR. INCAVG.EQ.3! ) XXLN=OMOVE 
CALL ASC ALE (TEMPX, XXLN, 2,1 , 10. 0) 

CALL aSCALEITcMPY, YLN, 2, 1, 10.0) 

X (NP* II =0. 

X(NP*2 J=TEMPX(4) 

XSC = T EMP X ( 4 I 
Y I NP* 1 J = 0. 

YlNP*2)=TcMPYI4) 

YSC=TEMPY (4» 


drag the x+y axes. 

IF( IUN1TS.NE.0) GO TO 69 5 

1F( ITEST .EQ. 0 .ANU. YMUIT .EG. .0011 CALL A XE S ( 0. , 0 . ,90 . , YLN , 0 . , 
1 YSC , l . * 0 . , 10HS TRESS, XSI, .25, 10) 

1F(ITEST .tQ. 0 .AND. YMULT .EG. 1.) CALL A XE S I 0. , 0 • , 90. » YLN , 0. » 
l YSC , 1 . ,0. ,10HSTRESS,PSI, .25, 10 J 

I F ( ITEST .EQ. 1 .AND. YMULT .EQ. .OOL) CALL AXES ( 0. , 0 . , 90 . , YLN , 0. , 
1YSC , 1 . ,0. , 16HSHEAR ST RES S , KS I , .25 , 16 1 

I F( IT EST .EQ. 1 .AND. YMULT .EC. L.) CALL A XE SI 0. , 0 . , 90. , YLN , 0. , 
1YSC, 1. , 0 . , 16HSHEAR STRES S , PS I , .25 , 16 I 

l F ( I 1 NDA V G «NE .2 J .AND. (ITEST .NE. 1 .OR. INOAVG .NE. I))GO TO 68 

IF( ( ITEST. EQ.O .OR. INOAVG. EC. 0 .OR, INOAVG. EQ, 3) .ANO. XMULT. EQ.l 
1) CALL AXESIO. , 0. , J. , XLN, 0., XSC, 1.0,0., 7HSTRAI N ,.25,-121 

1F1 ( ITEST. EQ. I .AND. ( INOAVG. EQ.l .OR. INOAVG. EQ.2J ) .AND. XMULT. EQ. 
Ill CALL AXcS(0.,0.,0.,XLN,0.,XSC, l. , o . , 13FISHE AR STRAIN ,.25, 

2-18! 

i F ( XMULT .EQ. 11 GO TO 69 

CALL AXESIO. ,0.,0., XLN» 0. »XSC ,-i. ,0. , 1H » .25,-11 
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i Ft ITEST.EQ.O « OR- l INOAVG.EQ-.O.OR. INDAVG.EQ.3) ICALL NOTATEC XLN/2. 
l,-.7,.25,18HSTRAIN,MICRO IN/IN.O. .18) 

IFC ITEST.EQ.l .AND. (1NOAVG.EQ.1 .OR. INOAVG.EQ.2 ) I CALL NOTATEC 

1XLN/2. ,-«7».25»24H5HEAR STRAIN, MICRO IN/IN, 0., 24) 

GO TO 69 

68 CALL DASHLNtO. ,0., XLN. 0. ,XLN) 

PLACc= IXSIZE/4 

CALL DASHLNCPLACE ,-.6, PLACE ,0.,.6) 

CALL DASHLNC PL ACE + I . 6 , PL AC E +1..0...6) 

CALL PaRR0W(PLACE-.3,-.30, PLACE, -.30, 1) 

CALL PAR ROW (PLACE+1.3,-.30,PLACE+l.»“.3Q» 1) 

CALL DUMBER I PLACE*. 2, -.3 7, .14, XSC/XMULT ,0. ,4) 

I F ( I TEST .EQ. 0) CALL NOTATE ( XLN/ 2. , - . 8 ,. 25 , 7HSTRAIN ,0.,12) 

IFUTcSI .EQ. 1) CALL NOT ATE C XLN/ 2. .8 ,. 2 5, 1 3HSHEAR STRAIN , 

10., 16) 

GO TO 69 

695 1F( ITcST.fcW.O .AND. YMULT.EC.l.) CALL AXES ( 0. , 0. , 90. , YLN, 0. , YSC , 1. 

1.0. , 9HSTRcSS, PA..25, 9) 

IFUTtST.EQ.J .AND. YMULT.EQ. .001) CALL AX ES ( 0 . , 0. , 90. , YLN , 0. , YSC 
1 , L. , 0., 10HSTKESS.XPA..25, 10) 

IF< ITtST.EQ.O .AND. YMULT.EQ. L0.**C-6)J CALL AXE S( 0 . , 0. , 90 . , YL N, 0 
1. ,YSC, 1. ,0. ,10HSTRESS,MPA , .25, 10) 

IFC ITlsT. EQ.O .AND. YMULT.EQ. 10. **(-9)) CALL AXE SC 0 . , 0. , 90. , YLN , 0 

1. , YSC,i. ,0.,10HSTRESS,GPA,.25,10> 

I Ft ITcST .cQ.O ,.*,N0. YMULT.EQ. 10. **l- 12 ) ) CALL AX ES ( 0 . ,0. , 90. , YLN , 

10. , YSC, 1. , 0. , 10H STRESS ,TPA ,.25,10) 

IFC ITLST.EU.l .AND. YMULT.EQ. 1) CALL AXES ( 0 . , 0. , 90. , YLN , 0. , YSC , l. , 

10. . 15,SHt AK STRESS, PA, .25, 15) 

IFC ITtST.EQ.l .ANO. YMULT.EC. .001) CALL AXES ( 0. , 0. , 90. , YL N , 0. , YSC 

1 . 1. . 0. . 16,6HEAK STRESS »KPA ,.25, 16 ) 

IFC ITcST. EU.l .ANO. YMULT.EQ. 10. **(-6)) CALL AXE SC 0. , 0. , 90. , YLN , 0 

1., YSC ,1. , 0. , 16 H SHEAR S TR £ S S , MP A , . 25 , 1 6 ) 

IFC 1 TEST . EQ. 1 .AND. YMULT.EQ. 10. **(-9)) CALL AXE S ( 0 . ,0. , 90 . , YLN , 0 

1., YSC,l. ,0., 16HSHEAR STRES S , GP A , . 25 , 16 ) 

IFC ITEST.EQ.i .ANU. YMULT.EQ. 10. **(-12)) CALL AX ES C 0. , 0 . , 90. , YLN , 

10., YSC, 1 • ,0. , 16HSHcAR STRESS,TPA,.25,16) 

IFC C 1N0A Vti.NE.2) .AND. CITEST.NE.l .OR. I NCAVG.NE . L ) ) GO TO 685 
IFC (ITcST. EQ.O .OR. INDAVG.EQ.O .OR. 1N0AVG.EQ.3) .AND. XMULT.EG.I 

I) CALL AXESCO. , 0. , 0. , XLN, O..XSC, 1.0,0., 7HSTRAIN ,.25,-10) 

IFC ( I TEST.EQ.i .ANu. ( 1NDAVG.E0. 1 .OR. INOAVG.EQ.2) ) .ANU. XMULT.EQ. 

II) CALL AXESCO. ,0.,0., XLN, C. , XSC , 1 . , 0 . 1 3 HS HE AR STRAIN ,.25,-16 

2 ) 

IFIXMULT .CQ.l )GU TD 69 

CALL AXfc SCO. ,0. ,0. ,XLN,0. ,XSC,-l. ,0. , 1H , .25,-1) 

IFC IT tST .EQ.O .UR. C I NOAVG .EQ.O. OR . [ NDA VG. E Q. 3 ) )C ALL NOTA TE C XLN/2 . 
l,-»7,. 25,1 8H STRAIN, MICRO P/M ,0.,I8) 

IFC ITEST.tQ.l .AND. (INOAVG.EQ.l .UR. INDAVG.EQ.2 ) ) CALL NOTATEC 
lXLN/2.»-.7,.25»24H$HEAR STRAIN, MICRO M/M ,0.,24> 

GU TO 69 

685 CALL OASHLNCO. ,0. ,XLN, 0. ,XLN) 

PLAC£=lXSlZE/4 

C„LL DaShLN(PLACE,-.6, PLACE,0.,.6I 
CALL UASHLNC PLACE +1. , - .6 , PLACE+1 . ,0 . , .6 ) 

CALL PaRROWCPL ACE-.3,- .3,PLACE,-.3, 1) 

CALL PARROWCPL AC E+ 1.3, -.3, FLACE+1. ,-.3,1) 

CALL NUM8ERCPLACE + .2 ,-.37, . 14, XSC / XMUL T , 0. ,4) 

IFC I TcST.EU.O) CALL NO TA TE ( XLN/2 . 8 , . 25 , 7HSTRAIN ,0.,10) 

IF( ITLsT.EQ.l) CALL NOTA TE ( XLN/2. , 8 , . 25 , 1 3HSHEAR STRAIN ,0., 

1 lo) 

69 CALL UGRIUCO. ,0. ,1., 1. , I XS I2E , 1YS I ZE ,0. 1 ) 

ROUNT = ROUNT + L 



APPENDIX B 


C 

C LOAD X ARRAY, PRINT, AND PLOT. 

1F< ITEST .E Q. 0 .AND. INDAVG .NE. 2JCALL COMP I ND I IG AGES, HOD ,NRLSLP 
l,NRUSLP»JAC,NRL»NRU»NP,X»Y, XMULT , I NDA VG , I XSIZE»DTA*YLN* ISPEC»YMULT 
2, 1 UN ITS) 

IF! ITEST ,EQ. 0 .AND. INDAVG .60. 2)CALL COMP AVG ( I GAGES, MOD, NRL SLP 
l , NRUSLP, JAC,NRL,NRU, NP , X, Y , XMULT , DTA , YLN, 1 SPEC ,YMULT , I UNI TS I 
IFIISPtC .tQ. 1 .AND. ITEST .EG. L .ANC. UNOAVG .EQ. 0 .OR. 

1 INDAVG .EQ. 3 ) (CALL TOR I NO 1 1 GAGES , MOD , NRL SI P, NRUSLP , JAC , NROS , NRL , 
2NRU,NP,X,Y, XMULT, INDAVG, IX SI ZE,DTA,YLN, XSC, RMAXSN ,YMULT, I UNI TSi 
IFIISPeC .EQ. I .AND. ITEST .EQ. 1 .AND. (INDAVG .EQ. I .OR. 

I INDAVG .EW. 2) ) CALL TORAVG ( I GAGES , MOD , NRL SLP, NRUSLP, JAC , NROS , NRL, 
2NRU,NP,X,Y, XMULT, INDAVG, I XSI 1 £ , DT A , YL N , XSC , YSC, FAC, YMULT, AL6NGTH , 
3RLc NG TM , I UNIT'S) 

IF( ITEST. EQ. I .AND. INDAVG. EQ. 2) KOUNT = KOUNT+ 1 
I F( MUD. cQ.2 ) KOUNT = KOUNT ♦ l 

C 

C GO Tu NEW GRAPH, NEW RUN. 

50 IFIMURE.EG.O .AND. MD.EQ.9999) GO TO 1000 
CALL CmLPLTULN+6. ,0. ,-3) 

PRINT 1200 

1200 FURMATUHl) 

IF ( MuKL. EQ.O) GO TO 55 
KNTRUN=K.NTkUN+ l 

IFl ITEST. EU. U) RcaDI 5, 700) NRL ,NRU, INDAVG, I GAGES , MOD , NRL SL P , NRUSLP , 
1NUPLT ,MQKt 

700 FORMAT l 3 13,. 412,313,101 2, 17X, 1 1) 

1 F ( ITtST.EQ.l) READ(5, 701) NRL, NRU, NROS , 1 ND AVG, MOD, NRL SLP, NRUSLP, 

1 NOPLT .MORE 

701 FORMAT! 713, 1312, 22X, 11) 

NP = NR U-NRL + 1 

GO I J 32 

1000 CALL CALPLTIO. ,0. ,9V9) 

PRINT 1025, KOUNT 

1025 FORMaK/// LH * NUMBER OF FLCTS = * 13) 

PRINT 1026, I XS 1 2c , 1 YS IZE 

1026 FORMA T ( 1 H *l£NGTH OF X AXIS = * 14 * INCHES* / 1H ^LENGTH OF Y AX 
US = * 14 * INCHES*) 

END 


SUBROUTINE RDTAPEl JAC , NR L , N RU , MD , DT A ) 

C 

C READS BECKMAN DATA FROM MAGNETIC TAPE IN FORMAT 2. 

C 

DIMENSION AC 2 1 99 J , DT A ( 4 1 5 , 20 ) 

PRINT 198 

198 FCRMmTUHO,* FRAME LOAD SG-l SG-2 *) 

N=2 

NP=NkU-NRL*l 

100 KEADt 30) MO, IF M, ISP, I FRA, TIM, TV, I FAC, I T ST , MRN , MSN , DC6 , DC 7 , 0C8 , DC9 , 
1UU0, A^2 
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c 

C CHECK FOR END OF RUN. 

1F(MD .tQ. 8883 .UR. MC .EC. 9999) GO TO 199 

PRluT EVERY FIFTH GR TENTH FRAME. 

IHINP .Lt. 100 .A NO . N/ 5*5 .EQ. N) .OR. (N/10*1G .EQ. N)) PRINT 
1 10 1 • IFRA, lAC2(i\> ,K=1,3) 

101 FORMATI1HO, l l2,3Elc.4) 

STORE DaT« IN AkRAY DTA. 

LKJ 10 0 J = 1 * JAC 
103 OTA ( iFKAt 0) = AC2( J) 

N = N + 1 
GO TO 100 
199 RETURN 
End 


SUBROUTINE RGCELLI JAC ,NRL ,NRU,MD, DTA) 

READS UATA FROM uATA FILE CN OATA CELL — OATA ON TAPE30 

DIMENSION UTA(415,20) 

PRINT 193 

l 9a FORMAT ( IH0»* FRAME LOAD SG-1 SG-2 *) 

N = 2 

NP=NRU-NRL+l 

A 00 RtAUl 30, 401) I T S T , MRN , MD , I FR A , I S Q 

401 FuRMATmi2,24X,Il2,l8) 

C CHlCK FOR END f)F RUN. 

IFIMU .EQ. 8888 .OR. MO .EQ. 9999) GO TO 410 
C READ OATA INTO AKRAY OTA. 

DO 411 J = l, 18, 6 
J J=J+5 

411 KcAUI 30 , 402 ) ( OTA ( J f RA , K ) , K= J, JJ ) , I SQ 

402 F0RMATI6E 12.4, 18) 

C 

C PRINT EVERY FIFTH UR TENTH FRAME. 

IFKNP ,LE. 100 .A NO. N/5*5 .EQ. N) .OR. (N/10*10 .EQ. N)) PRINT 
1403, IFRA, IOTA! IFRA, K) ,K=l, 3) 

403 FORMAT! IHO, 112,3012. 4) 

N=N + l 

GO TO 400 
410 RETURN 
END 
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SUBRJUT I NE LOAOY (THICK ,WIDTH»SCFAC, I SPEC* I TEST , NRL, NRU,Y * DTA* FAC* 

1 YMIJIT , I UNIT S ,KNTRUN) 

C COMPUTES STRESSES AND LOADS Y ARRAY FOR EACH SPEC * TEST TYPE. 

c 

DIMENSION YI400) , DTA ( 4 15 . 2C I 
IF(KNTkUN.NE.O) GO TO 19 
Pl=3. 141S9 

[FUSPaC .EQ. O) fac=thick*wioth 

I F ( I SPEC .EQ. I .AND. ITEST .EQ. 0) F AC=P I * ( W I DTH-THl CKJ * THIC K 
1 F l I S PEC .63. 1 .AND. ITEST .EC. I) F AC=2 *P l * ( ( {WIDTH- THlCK)/2) **2 
1J*THICK 

1 F I ISP6C.tO.OJ PRINT 191, THICK, WIDTH 

191 FURFUTI IHu,* THICKNESS = * , £ 12 .4 , 6X , * WIDTH = *,E12.4) . 
IFIISPEC.6U.il PRINT 192, THICK, WIDTH 

192 FORMAT ( 1 H 0 , * THICKNESS = *,E12.4,6X,* OUTER DIAMETER = *,E12.4) 

C IF STRcSSlS TOO LARGc, DIVIDE BY 1000. 

L9 YMULT=l.O 

RMmXL D-ABSIDTA (NRU, l ) ) 

00 196 I =NRL , NRU 

196 IFIaBaIDTAI I , l ) ) .GT.RMAXLD) RM4XLD=AB S ( DTA ( 1,1) > 

RMAXSS=KMAXL3/FAC 

IMRMAXSS .GT. 1000) YMUtT = .001 

1 F ( 1 UNITS. EQ.O I GU TlJ 19 3 

I F ( R MAX S S* . 00 1 .GT. 1000) YMUL T= l 0 . ** ( - 6 ) 

IFIRMAXSS*(10.**(-6) ) .GT . 1000) YMULT= 10 . **(-9) 

IFIRMAXSS*! 10. ’**1-9) ) .GT. 1000) YMULT= 10. -12) 

C PUT STkcSSES IN Y ARRAY. 

193 DO 200 1=NKL,NRU 

2 00 Y< l-NRL + 1 IMDTAI 1 , 1 ) / F AC ) * YM OL T*3CF AC 
RETURN 
END 


FUNCTION POLY E 1 IX»M,C) 
UATA BIG/0' 1 77777777777/ 
DIMENSION CIM) 

I F ( M- 1)10*11,12 

12 N = M-1 

F UL Y 1 1 — C l 1) 

00201=1 »N 

20 POLYtl=X*POLYEl+C( I «-l ) 
RETURN 

ID POLYtl=oIG 
RETURN 

11 POLYtl-CI 1) 

RETURN 

END 
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SUBROUTINE R DC ARC l JAC , NR L , NRU»MD» DT A ) 

REAUS DATA FRuM PUNCHED CARDS. 

DIMENSION UT A I 41 5 , 20 J 
PRINT 1 9 d 

198 FORMAT ( 1H0,* FRAME LOAD SG-1 SG-2 *) 

N=2 

NP=NKU-NRL+l 

500 DO 511 0 = 1 t J AC ,6 
J J - J +■ 5 

IFIJ.fcQ. U READ! 5, 501) If RA , l CT A 1 1 FRA, K) , K= J , J J 1 » MO 

501 FURMAT (I 4»6E12.<*, 141 
IF(J.Nt.l) RLAD!5,502) I DT A ( I FRA , K I , K= J , U 0 ) , MD 

502 FUKMATI6E12.4, 14) 

511 CONTINUE 

C CHECK FOR END UF RUN. 

I F ( M 0 .HO. 8388 .OR. MD .CO. 9999 1 GO TC 510 

C PRINT liVERY FIFTH OR TENTH FRAME. 

I F ( ( N P .LE. 100 .AND. N/5*5 .EU. N) .OR. <N/L0*10 .EO. N) ) PRINT 
1503, IFRA,(DTA( I FR A , K ) ,K= 1 ,3 ) 

503 FORMAT! I 12,32 12.4) 

N=N+ 1 

GO TO 50 0 
510 RETURN 
END 


SJBROUTINc CUMPIND! I GAGES, MCC.NRLSLP, NRUSLP , J AC , NRL , NRU ,NP , X , Y, 

IX MULT , INDAVG, I XS IZ t, UT A , YLN , I SPEC , Y MULT , 1 UN I TS ) 

C 

C LOADS X ARRAY FUR COMP TEST, INDIV GAGE PLOT, AND CALLS PRTPLT. 

C 

01 MANSION I GAGES (4) ,X (400) , Y< AOOI ,DTA( 415 , 20) 

CoMrtJN UMOVE , I PLOTS , NOPL T ( Id ) , XLN , A Y , BX , I UT S 
DO 122 I G AG=2 , J AC 
DU 124 K = l, 13 

124 IFl IGAG-1 .tO.NOPLT (K) ) GC TC 122 
Do 121 I = NRL, NRU 

i 2 1 XI 1-hRl+ l ) =AiJS ( OTA 1 1 , I GAG) ♦XMULT I 

CALL PRTPLT ( I SPEC, 0, INOAVG , I GAGES, I GAG , NR L , NRU , NP , X , Y , YLN > 
l F ( I j AG .LT. JAC) CALL C A L P LT l CMOV E , 0 . , -3 ) 

122 CONTINUE 

IF ( INOAVG .EU. U) GO TC 12J 
DO 12 5 K = 1 , 18 

IFIJAC-I .cU.NUPLT(K) ) GO TC 126 

125 CONTINUE 

CALl CAL PL T ! OMQVt , 0. 3 ) 

126 CALL CuMPAVG < l 0 AG E S , MOC , NRl SLP , NRU SLP , J AC , NRL . NRU , NP , X , Y , XMULT , OTA 
1 , YlN, l SP EC , Y MUL T , I uNI T S) 

123 RETURN 
END 
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SUBROUTINE COMPAVGIIGAGES, RCO , NRL SLP , NRUSLP , J AC , NRL , NRU, NP , X , Y, 

1 X MULT , OT A , YLN, ISPEC.YMULT. 1UNITS) 

LOADS X ARRAY FOR COMP TEST, AVG PLOT, AND CALLS PRTPLT. 

IF MOO= 1 , CALLS SLP. 

DIMENSION IGAGES 14) ,X(4GO) , Y ( 400) , DTA ( 4 15 * 20) 

COMMON UMOVE , I PL CT S , NO P L T ( 1 8 ) , XL N , AY , BX , I UT S 
FIND NUMBER OF GAGES TO BE AVERAGED. 

N=0 

DO 101 I =1,4 

IF(IGAGESU) .NE. 0) N = N«-1 

131 CONTINUE 

FIND AVERAGES *■ PLOT. 

DO 132 1=1, NP 

132 X ( l ) = 0 

DO 135 J=1 ,4 

IFIIGAGES(J) -EU. 0) GC TC 135 

00 134 I = NRL , NRU 

I GAGE S J= I GAGES I J ) 

134 X ( I -NRL + 1 ) =X l l-NRL + lUABSICTAl I , i GAGESJ +1 ) * XMULT/N) 

135 CUNTINUE 

CALL PRTPLTI l SPEC, 0,2, IG AGES, 0, NRL , NRU, NP,X,Y,YLN I 
IF MuD=l, FIND ♦ PRINT SLOPE. 

1 FI MOU .tO. 0) GO TO 139 

IFlMOU.cU.2) CALL PROP L l M ( X, Y , NRL , NRU »NP » XMULT » YMULT , I UNI T S ) 
IFlMUu.tQ.DGJ TO 139 
DO 1331 I =’NRL SLP ,NKUSLP 
X(1-NRLSLP*U=XC l-NRL+1) 

1531 Y( I-NRLSLP + 1 ) = Y ( I-NKL+1) 

NPSLP=NRUSLP-NRlSLP+1 

CALL SLP (X, Y,XBR,Y6R,EX,NPSLP,XMULT, YMULT ) 

PRINT 1*8, EX, NRL SLP.NRUSLP 

136 FURMATl// 5X , * t = * » fc 12 . 5 ,4X ,* FOR FRAMES 1 *, 14,* THR0UGH*»I4 //) 
139 RLTURN 

END 


SUBROUTINE TOR I NO UG AGES ,MCD, NRL SLP, NRU SLP, JAC, NRCS.NRL, NRU,NP,X,Y 
l , XMUL T , I NDAVG, I X S U b , OTA , YLN , X SC , R MAX SN , YMULT , I UN I TS I 

LOADS X ,,RRAY FOR TUBc SPEC, TORSION TEST, INDIV GAGE PLOT, AND CALLS 
PRTPLT . 


160 
1 b6 

161 


DIMENSION IG AGES I 4), XI 400) , Y 14001 ,0TA(4 15,20) 
COMMON OMuVt , I PLOTS, NL'PLT ( 16 ) , XLN , AY, BX , IUT S 
l G aG= 3 

DO 166 *. = 1,16 

IFl IG aG-1 .EO.NOPLT IK ) ) GC TO 167 


DU 1 6 1 i ~\kL , NRU 

X ( I -NK L+U = ARS(CTA( l , l GAG ) *XMULT ) 

CALL PRTPLTI 1, 1, I NDAVG, I GAGES, I GAG , NRL , NKU, NP , X , Y , Y L N ) 
1 R ( I GaG/ 3*3 .cQ. IGAG) GO TC 162 





O DO J 
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I Ga b = IGAG + 2 

IF( IGA6.GT.3*NRGS+i.) G C TO 164 
GU TG lt>3 

162 l GAG- I GA G“ i 

163 L*LL CaLPU < LWOVt, 0. t-3) 

IR lGnb,Nc.-J « aNC • 1 GAG/ 3*3 . cQ • IG AG ) 

GU TU 16 0 

167 IF ( 16 AG/ 3*3 • EO • 1UAG) GC TG 166 
IGAb- lb A G + 2 

IK IGaG.GT. 3*NROS+ i) GC TC 164 
GO TO 160 

168 I GAG= 1GAG - 1 
GU TO 16 0 

164 IK I NDA VG • E 0 « 0) GO TC 16 5 
DO 1641 K= 1 1 18 

If! 3*NROSK.Eg .NCPLTIK) ) GC TG 1642 


CALL CAL PL T ( l.,0.,-3) 


1 o4 l CONTINUE 

CALL CaLPLT(UM0Vc, 0.,-3I 
l 6*t2 DO 166 I = NRL , NRU 

165 X ! 1-NRL* 1 1 = l ABS! DTAl 1,2) > + ABS( DTA l 1,4) ) )*XMULT 

CALL PRTPLT(1,1*INCAVG,IGAGES»-1* NRL ,NPU,NP,X,Y,YLN) 

C IF MOD=l» FIND + PRINT SLOPE. 

I F ( MUU .LG. 0 ) GO TO 165 

IKMOD.EO.2) CALL PR OP L I M ! X, Y , NRL , NRU ,NP , XMUL T, YMULT , I UN l T S ) 
IF t hUD. fc Q.2) GU TC 165 
DO 1631 I =NRLS LP , NRUSLP 
X( 1-NRLSLPU) = X! I-NRL+1) 


1531 Yll-NRLSLP+U=Y( l-NRL + 1) 

NPSLP=NRUSLP-NRlSLP+l 

CALL SLP<X,Y,X8R,YBR,EX,NPSLP»XMULT,YMULT) 

PRINT 139,EX,NRLSLP, NRUSLP 

135 FORMAT!// 5X r * G = * , E 12 . 5 , 4X , * FOR F RAME S* , I 4 , * THROUGH*, 14 //) 
169 RETURN 
END 


SUBROUTINE TQRAVGl IG AGES , PCD , NRLSLP , NRUSl P , J AC , NROS , NRL , NRU, NP » X , Y 

1 , XMULT , 1.NDAVG, IXSIZE, DT A , YLN , XSC , Y SC , FAC , YMUL T , AL ENGTH, RLENGTH, 

2 I UN I T S ) 

LOADS X AkRAY FOR TUBE SPEC, TORSION TEST, AVG PLOT, AND CALLS 
PRTPLT. IF M0D=1 » CALLS SLP. 

DIMENSION !GAGcS(4) ,X(4001 ,Y!400J , DTA ( 4 15 , 201 , HOLDX ( 4 > ,HOLOY { 4) 
CUMMuN OMOVE, I PLOTS, NO PL T ! 18 ) , XLN , AY, BX , I UT S 
DU 171 1 = 1 , NP 

i7i xm=o 

N=0 

DU 172 0=1, NROS 
DO 161 K = l , 1 8 

IF! 3* J-2 .Et).NUPLT(K) .OR. 3* 0 . EQ .NUP LT ( K) ) N=N*1 
IF! 3*0-2.EU.NOPLT(X) .OR. 3* 0. 6Q.NOPLT! K1 > GO TO 172 
161 CONTINUE 
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lc 3 XU -NR L* l ) = X( I -NRL + U+ ( AfiS ( CTA( I, 3*J- 1 1 )*ABS < OTA( )*XMULT 

uz continue 

OU lo2 I = NR L * NRU 

162 XC1-NKL+ 1 >-X( I-NKL-*-! ) / (NF.OS-N ) 

CAll PR TP L T ( i i 1 » INCA VG » I GAGES, 0, NRL ( NRU f NP , X, Y , YLN) 

C IF MOU= 1 ? F I NO «■ PRINT SLOPE, 

lfCMUD .GQ. 0 1 GO TO 175 

IFIMUD.EQ.2) CALL PRGPL 1 M ( X, Y ,NRL , NRO fNP, XMOLT, YMUL T , 1 UN I TS 1 

IFIM0U.EQ.2JG0 TC 175 
DO 15 31 I =NRLSLP f NRUSLP 
XII-NRLSLP«-1)-X( l-NRL+l) 

1531 YU-NKLSLP+1 ) = Y( I-NRL+l) 
npslp=nruslp-nklslp+ l 

CnLL SLP (X, Y , XBM ,Y6R» EX,NPSLP»XMULT, YMULT ) 

PRINT 136, LXf NRlSLP, NRUSLP . _ • ... 

136 FORMAT!// 5X,* G = *,E12.5,4X,* FOR FRAMES#, 14,* THRG'JgH*»I4 //) 
C IF InUAVG = 1 » PUJT lWlST CURVES. 

175 IMINDAVG . Me . 2)GG TO 179 
P I = 3. 14 l 59 
M l N = J AC- 3 
MAX = M l N* 3 

C MAK.E Ntirt GRAPH. 

XLh=XLN 

IFIXlM.LT. 9. ) XLM=9. 


CALL CALPLTIXLM+6. ,0. »-3) 
HDLOY(l) =0. 

HQLOY ( 2 J =DTA ( NRU , 1 I 

CALL ASCALfelHOLUY,YLN,Z, l, 10.0) 

YSC2=H0L JY ( 4) 

Y ( NP + 2 ) = Y 3C2 


HOLDXl i) =0. 

HGlUX l 2 J = (ABS(OTA(NRU,'MIN) ) +• ABS ( DT A( NRU, MIN+11 ) -ABS COT A ( NRU, 

1M1N+2 U-A6SIUT A (NRU, MIN+3) ) )/ I 2*AL ENG TH*R LENGTH ) 

If U UNITS. EQ.O) HGLUX (2 ) = HCL0X(2) *160. /PI 
CALL ASC ALc ( HGLUX , XLN , 2, 1 , 10 « 0) 

X SC 2= HUL OX ( 4 ) 


XlNP+21 =XSC2 

IFl I UN ITS. EQ.O) CALL AXtS(0.,0.,90.,YLN,0.,YSC2,-l.,0.,l2HTORQUE,I 
LN-L 6, .25 , 12) 

1 F ( 1UN IT 5 . tQ . 1 ) CALL AX£SlO.«0.,90.,YLN,0«»YSC2i-l.»0.*L5HTQRQUE»N 


1 c *vT l) N—M , « 25 t 15) 

CALL AXt $< 0 . t 0. * 0# ?XLN 1 0« » XSC2*1« ? 0* t iH » 0*rl) 
CALL NUM BE R ( — l #5/28* »■ *25 » *1875, Q« i 0 • t — 1 ) 

DO 1771 1=1, IXS1ZE 


PLACE - I +0 • 

1771 CALL NUMBER (PLACE- 15./ 56 * , - .25, . 1 £75 , XSC2*PLACE,0.,4) 

IH IUUIT S.EQ.O) CALL NCTATEC XLN/2 . - U 3 1 -. 75 f .25 , 12HTWIST , DEG/ IN , 0 . 

1 I M IiiMl T $ *tQ. l ) CALL NCTATE(XLN/2.-L.3f-. 75, .25, llHTlrflST , RAD/M, 0. , 


ill) 

IVSIZE=YLN 

C At L DOR 1 L) ( 0« ,0* , 1 • , I • flXSIZE, IYSIZE,0. L ) 


C 


LOAD X ARRAY. 

DO 176 I =NKL , NRU 

X 1 1 -NRL * 1 ) = (A0SiCTA( I , MIN ) ) ♦ ABS < DTA ( 1 , M 1 N+ 1 ) ) -AB S l DT A ( I , M I N+ 

12)J-AdS( DTA( I,M1N + 3J > ) / < 2 . *A LENGTH*RL E NGTH ) 
l F ( IUN I T S • tO ♦ 0 ) X ( I — NR L+ 1 ) * XU-NRL + l) *180. /PI 


176 CONTINUE 



® poor 
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C LOAD Y ARRAY. 

DO 177 1 =NKL , NRU 
177 Y( I-NRL+l )=ABS< CTAUt 1 >) 

C PRINT + PLOT 

CALL PRTPLTU,l.INCAVG,lGAGES,-2, NRL , NRU, NP ,X , Y ♦ Y LN I 
1 F ( MOO .N c • 1 ) GO TO 179 
DO 1532 I =NRL$LP, NRUSLP 
Xtl-NRLSLP+l)=Xt I-NRl+1) 

1532 Y { 1— NRLSLPU l=YII-NRL+l) 

NPSLP=NRUSLP-NRLSLP*l 

CALL SLP(X,Y,XBR r YBR,tX,NPSLP»XMULTf YMULTJ 
PRINT 154 »EX,NRlSLP,NRLSLP 

159 F OR Mm T (//5X*SL0P£ OF CURVE = * E12. 5, 4X,*F0R FRAMES* 
1 14 * THROUGH* 14 //I 

179 RETURN 
END 


SUBROUT I ME PKTPLTl IS PEC, t T EST , I NDAVG , IGAGES , I GAG , NRL , NRU , NP, X,Y, 
1 YLN I 

PRINTS + PLOTS POINTS FROM X,Y ARRAYS. 

DIMtNSlUN 1 GAGES ( 4 I , X(40OI» YI400) 

DIMENSION N0PLT2 UB) , 1 FRAM1400) 

COMMON UMOVE , 1 PLOTS, NOPLT ( 10 ) »XLN, AY, BX, JUTS 
1XM = 1 
1 YM = 1 

IFRAM(l) = NRL 
DO 25 K = 2 , NP 

25 IFRAM(K) = IFRAM(K-l) +• 1 
IF( fcjX.cO . 11 1XM = 1 
I F I OX . cQ . I J**fa I 1 XM = 2 
I F( lUTS.cQ.il GC TO 20 

U.S. CUSTOMARY UMTS 

1 f ( AY .£(].! ) IYM = 1 

IFUY.cQ . .001 I IYM = 2 
G J TO 21 

C S I UN I T S 

C 

20 IFUY.EQ.l) IYM = 1 
IFlAY.cQ. .001) IYM = 2 
IFIAY.cQ. 1.E-J6) IYM = 3 
IFUY.EQ.l. £-09) I YM = 4 
1FUY .EG. l.c-12) IYM = 5 

21 CONTINUE 
RL=NKL*0. 

KU=NRU+0 . 

I F ( IGAG .LE. 01 GO TO 183 
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JGAG- IGAG-1 

PKINT 181 , JGAG, NRL, NRU 

181 FORMA T( 1H1 ,* STRESS-STRAIN POINTS FOR SG-*,i2»* FRAMES *»I3r* TH 
1 ROUGH * » 1 3 1 

CALL NuTATt < L. ,YLN*.35,. 1 , 3HSG-»0. , 3» 

GAG= 2 GAG- 1 • 

CALL NUMBER! 1. +1*8/7., YLN*. 35,. 1, GAG, 0. U 
GO TO 189 

183 IFUTEST .NE. 0) GO TO 186 
PRINT 184, 1 GAGES , NRL , NRU 

184 FORMAT! 1H1 ,*STR£SS-STRAIN POINTS FOR AVG OF GAGES *»4I3,* FRAMES 
i*,I 3,* THROUGH * , I 3) 

CALL NuTATc U. ,YLN*.35,. l,6HAVG 0F,0.,6) 

N=G 

00 L845 1=1,4 

lF(lGAGtS(ll .EQ. 0) GC TO 1845 
UM = I GAGES! I ) 

N=.N+ 1 

CmLL NUM BER ( 1 • F3»6/ 7 • +K* 1 .6/ 7. » YL N+« 35, . 1 , UM * 0. ,— 1) 

1845 CONTINUE 
GO TO 189 

186 IH1GAG . NL . -2>G0 TO 188 
PkINT 187, NRL ,NR J 

187 FURMATIIM,* TORQUE TWIST POINTS, FRAMES #,I3,* THROUGH *»I3) 

GU TO 189 

188 IF! IGAG.NL. OJ GG TO 1882 
PRINT 1881, NkL.NRU 

1681 FurMAT ! IH’1 ,* AVG STRE SS- STRA I N POINTS FOR FRAMES *,I3»* THROUGH *, 
113) 

X = l 

00 10 N= 1,18 

LF!nlSPLTIN).EU.O) GU TC 10 
NUPLT2(M=n0PLT!N) 

K.=K + 1 

10 CONTINUE 
K. =K— 1 

lFlK.EQ.O) GU TO 15 
PRINT 11 

11 F OR MA Kiri ,* u AGES OMITTED ARE *> 

PRINT 12, !NJPLT2(N1 ,N=l , K ) 

12 FORMAT! 1H* ,22X , 181 3J 
15 CONTINUE 

CALL NUT AT l! 1. , YLN + . 35 , . 1 , 8HAVG PLCT,0.,8) 

GO TO 189 

18 82 PRIM 18 83, NRL, NRU 

1883 FORMAT! IH1 ,* STRESS-STRAIN PCINTS FOR 5UM OF GAGES 1 ANU 3, FRAMES 
1 *, 13»* THROUGH ^,131 

CALL NUTATE! 1. ,YLN + . 35,. 1, UHSG-1 ♦ SG-3*0.,li) 

189 Continue 

IF! IG4G.LQ.-2) GO TO 216 

PRINT 900 


900 FORMAT! 1H0, IX, 4 ( 28HFR AME STRESS STRAIN )) 

IF! I UTS . EQ. 1 ) GO TO 22 
GO TO 1200,201) , IYM 


200 PRINT 901 

901 FGRMmUIH , 1 X , 4 ( 2 8H 

PSI 

) ) 

GO TU 203 



201 PRINT 902 

902 FURMATUH ,1X, 4 ( 2 8H 

203 GU TO 1204,205 ) • I XM 

KS I 

J ) 
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204 
9U3 

205 

904 

905 

22 

206 

906 

207 

907 

208 

908 


209 

909 

*10 

910 
212 

213 

911 

214 

912 

913 


216 

C 

L 


PRINT 903 
FORMA T ( l H ♦ » IX, 

4(23H 


IN/IN 

) > 

GU TO 215 
PRINT 904 
FORMAT ( 1 H + » IX » 

<t( 23H 


MICRO 

) ) 

PRINT 905 
F ORMA T ( 1 H flX, 

4<2dH 


IN/IN 

) ) 

GO TO 215 
GO TJ (206,207 
PRINT 906 
FORMATUH ,1X, 

,203,209,210) , 
4(2oH 

I YM 

PA 


I J 

GJ TO 212 
PRINT 907 
FORMAT ( l H , IX , 

4I23H 

K P A 


) ) 

GO TO 212 
PRINT 90 0 
FuKMmTUH ,1X, 

4 ( 2 3 H 

MPA 


) ) 

GO TO 212 
PRINT 9U 9 
FORMaTUH ,1X, 

4 ( 2 3H 

G P A 


) ) 

GO TO 212 
PRINT 910 
FORMaTUH ,IX, 

4 ( 2 3H 

TP A 


) I 


GO TO ( 2 1 3 1 2 14 ) i IXM 
PRINT 911 

FUKMATl 1H+, IX, 4(2 8H M/M )) 

GO TO 213 
PRINT 912 

FORMAT ( 1H+, IX, 4(28H MICRO >) 

PRINT 913 

FURMATl 1H f IX, 4< 2dH M/M ) ) 

GU TO 215 

CONTINUE 


TORQUE TWIST DATA 


PRINT 915 


915 

FORMAT I IHO, LX, 4I28HFRAME 

TORQUE 

TWIST 

) ) 

916 

I F ( IUIS.EQ.1) GO TO 217 
PRINT 916 

FGkMATUH ,1X 4( 28H 

IN-LBS 

DcG/lN 

) ) 

217 

917 

GO TO 215 
PRINT 917 

FORMATUH ,1X, 4I28H 

N-M 

RAC/M 

) ) 


215 CONTINUE 

PRINT 191, M IFRAMIN) ,Y(N) f XINH tN =1,NPI 
191 FORMAT ( 1H ,1X, 4U4,2E12.4)> 

CALL NUTATE: C 1. t YLN+.2, *1 tl IHFR - ,0.tll) 

CALL NUMBER! l#+l«8/7«,YLN+*2»*ltRL?G*f— 1) 

CALL NUMBER! l#+4.d/7« , YLN* .2 t ♦ It RU t 0. 1 ) 

CALL LINPLKX, YfNPtltO rOt ItOI 

RETURN 

END 
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SUBRuUT 1NE SLP t X, Y , X BR , Y BR , S P t NU, X MU L T , YMSJL T1 
L 

C KUO T I <\lc TU DETERMINE LINE CF BEST FIT 

C Y- YBR = sP*( X - XBR 1 

C 

DIMENSION XUUO) ,Y(400I 
DU 20 1 = 1, MU 
X ( I J = X ( IJ /XMJLT 
20 Y ( I ) = Y ( I i/YMJLT 
XBR = X ( 1 > /NO 
Y BR = Y( 1 I /NU 
Du 10 JK =2, NU 
XBR = XBR *■ X( JK 1/Nu 
YbR = YBR ♦ Yl JKJ/NU 

10 CONTINUE 
C 

AN = XI 1 )* YU 1 
DO 11 JK = 2 i NU 
AN = AN + XIJKJ* Y(JM 

1 1 CONTINUE 

Au = xm* xm 

DO 12 JK =2 » NU 

AD = AO * X ( JK ) * X(JK) 

12 CONTINUE 

SR = < AN - NU* XBR* YBR 1 / (AC - NU* XBR * XBR! 

RETURN 

END 


SUbROUT 1 Ne PROPLI M( X , Y ,NRL ,NRU,NP , X MULT ,YMULT, l UNITS ) 
DIMENSION X14U01 , Y<4UG ) , ETAM4001 

DIMENSION YPS1 400! , YPTS< 40C, 11 ,W( 4001 , RES 10(400* 1 > »SUM< 1 1 , 

1 A( j, J1 ,b( 3, 1 1 ,C( 400, 1) j XP ( 400 ) 

COMMON UMCV£,IPLOTS,NuPLT< 18 1 ,XLN, AY, BX, I UTS 
DATAw/400*1.0/ 

OU 10 l = 1 » NP 
X( 1 l=Xl I ) /XMULT 
10 Y(1 1 = YI I 1/YMJLT 
00 150 J = 1 , N P 
YPTS( Jill = Y< J! 

XP(J1 = X(J1 
150 CONTINUE 

CALLLSwPOLiXP, YPTS ,* ,RLS I C , NP , SUM , 1 , A , b ,3 , C ,400, 3 1 
DO 151 J = 1 , NP 

151 ETAN(J) = b(c,il ♦ 2.*B( 3* 1)*XU) 

DO 152 KT = l f NP 

152 YPSIKT! = b { 1 , 11 * ai2,U*XCKT) ♦ B(3,ll *X<KT!*XIKT1 
PRINT 60 

60 FOKMmTUHI,* STRESS STRAIN TANGENT MODULUS*! 

PRINT 70, mm, X ( I I, ETAN( 1 1 ) , I =l,NP) 

70 FORMAT 1 1H ,31 g 12*4, 3X) ) 

CALL CALPlTIXLN*6. , 0. ,-3 1 

CALL BSC ALE (ET AN, 9., NP, 1,0., -1,0. ) 
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CALL BSC AlEI Yf 9.,NP»ltO.»-L«O.I 
[p{ IUN1TS .=Q.O I CALL AXES(0.»0.» 
1 » 19HT ANGENT MJUULU S » P S I » • 2 5 » 1*9 ) 
i F ( IUNiTS.NE.J 1 CALL AXESIO.fO.f 
1, IbPtTANGENT M0 UULUSfPAf.25f18) 

I p ( [UNITS .tD .0 1 CALL AXESIO.fO.f 
1ESS.PSI , .25,-101 
IF( I UNI TS.NE • 0 1 CALL AXESIO.fO.f 
1SS, PA, .25,-91 

CALL DGRIDIU.fO.fI., 1. ,9,9,0.11 
CALL LINPLT lY,fcTAN,NP , 1»0»C» 1*01 

return 

END 


90. ,9. , ETAN(NP*il ,ETANINP*21 , 1 . ,0. 

90*»9.»ETAN(NP+l) ,ETANlNP+21 , 1 « , 0. 
0.f9.,YINP+11 ,Y(NP+2) ,1.,0. , 10HSTR 
0.,9.fY(NP*1)fY<NP* 2),1. i0.i 9HSTRE 


SUbROUT INE LSQPOLlXtY,W,RESID»NiSUM»L»A ,b »M, C t NMAX ,MMAX1 
********* DOCUMENT DATE 08- 01 —b 8 SUBRUlTlNfc REV Is ED ub-01-oU ********* 

” UlMtNSIuN XINMAX1 ,Y (NMAX, L> ,RES IDlNMAA, L 1 »AlMMAX» MMAX1 , 

2blMMAA, L 1 » C l N M A X . M ) , SUMIL) , WlNMAX) 

C 

10 DJ 20 I = 1 , N 

2U III #11=1.0 

SO JO SO J=2.M 

9 0 DU SO I =1 # N 

SO Cl I F J) = C< I . J- 1 ) *X< I ) 

t»u ou loO i = 1 » M 

70 DU 100 0=1,* 

du AI I *01=0.0 

90 DO 100 R = 1 * N 

100 AU,J)=A( I fJ> *C(K, I >*C(K,Jl*W(lO 

10 5 JO ISO J = 1 1 L 

110 OU ISO 1=1 #M 

ItO til 1 . J)=U.0 

1J0 00 Iso K= 1 * N 

150 b ( 1 t j)=b( I #J) +C l K , 1 )*YlK,J)*w(K) 

CALL MATINV ( A# M # B# L# DET tRM , RE S 10 # u # UMAX , 1 S», ALE 1 
lbO 00 205 J = 1,L 
18i< SUflld)=U.O 
RK=M 

192 OO 19S K=1,M 

ClK,l)=Bl KK, J 1 
19S KR=KR-1 
198 UU 20S 1=1 fN 

RES Id I , J )=P0LYF1 IX t I 1 , M , Cl-Y { I ,0 J 
20s SUM l J»=SUM(JJ*RESIDll , J ) **2*W I I 1 
210 «.cTURN 
END 
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DESCRIPTION OF FORTRAN VARIABLES 


The following list contains a description of the significant FORTRAN variables 
appearing in the program. The dimensions for each array are given by using the nota- 
tion A(m,n). The variables which are input are not presented here since they are 
described in the section "Input.’' 


FORTRAN variable 
AC2(99) 

B(3,l) 

DTA(415,20) 

E TAN (400) 

EX 

FAC 

IPLOTS 

KNTRUN 

KOUNT 

NP 

NPSLP 

PI 

RMAXSS 

X(400) 

XMULT 


XP(400) 


Description 

temporary storage for input data 

stores coefficients of the polynomial fitted to the data points 
stores input data for all channels 
stores tangent moduli at each data point 

slope of the line of best fit to a specified section of the curve 

number by which load is divided to get stress (refer to table I) 

number of curves to be plotted for this submission; determined 
by the computer from information given by the user 

number of plots that have been made for this submission 

number of graphs that have been made. This is printed at the 
end of output. 

number of points to be plotted 

number of points to be used in finding the slope of the line of 
best fit 

PI = 3.14159 

maximum stress to be plotted 
stores strain values to be plotted 

strains, read in units of in./in. are multiplied by XMULT; if 
the maximum strain is less than 0.1 in./in., XMULT = 10 6 
and strains are changed to micro in./in.; otherwise, 

XMULT = 1. (When data are expressed in SI Units, strain 
is in m/m and micro m/m.) 

temporary storage area used in computing the second-order 
polynomial which provides the least-squares fit to curve 
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FORTRAN variable 
XSC 
XSC2 
XXLN 
Y(400) 

YMULT 


YPTS(400,1) 

YSC 

YSC2 


Description 

units per inch on X-axis for stress -strain curves 
units per inch on X-axis for torque -twist curves 
distance over which strain values are scaled 
stores stress values to be plotted 

stresses in psi are multiplied by YMULT; if the maximum 
stress is greater than 1000 psi, YMULT = 0.001 and stresses 
are changed to ksi; otherwise, YMULT = 1. (When data are 
expressed in SI Units, stress is in Pa and MPa.) 

temporary storage area used in computing the second-order 
polynomial which provides the least-squares fit to curve 

units per inch on Y-axis for stress -strain curves 

units per inch on Y-axis for torque -twist curves 
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LANGLEY LIBRARY SUBROUTINE MATINV 

Language : FORTRAN 

Purpose: MATINV solves the matrix equation AX - B, where A is a square coefficient matrix and B is a 
matrix of constant vectors. The solution to a set of simultaneous equations, the matrix inverse, and the 
determinant may be obtained. 

Use : CALL MATINV (A ,N,B ,M,DE TERM,IPIV OT , INDEX, NMAX , ISC ALE) 

A A two-dimensional array of the coefficients. On return to the calling program, A~^ is 

stored in A 

N The order of A, 1 = N = NMAX 

B A two-dimensional array of the constant vectors B. On return to the calling program, X 

is stored in B 

M The number of column vectors in B. The expression M = 0 signals that the subroutine 

is used solely for inversion; however, in the CALL statement an entry corresponding to 
B must still be present 

DETERM Gives the value of the determinant by the formula DET(A) = (l0 100 ) ISC ALE (DETERM) 
IPIVOT A one-dimensional array of temporary storage used by the routine 

INDEX A two-dimensional array of temporary storage used by the routine 

NMAX The maximum order of A as stated in the DIMENSION -statement of the calling program 

ISCALE A scale factor computed by the subroutine to keep the results of computation within the 

floating-point word size of the computer 

Restrictions: Arrays A, B, IPIVOT, and INDEX have variable dimensions in the subroutine. The maxi- 
mum size of these arrays must be specified in a DIMENSION statement of the calling program as 
A( NM AX , NMAX) , B(NMAX,M) , IPIVOT(NMAX), and INDEX(NMAX,2). The original matrices A and B 
are destroyed. They must be saved by the user if there is further need for them. The determinant is 
set to zero for a singular matrix. 

Method: Jordan's method is used to reduce a matrix A to the identity matrix I through a- succession of 
elementary transformations l n , l n _ i, . . l\. A = I. If these transformations are simultaneously 

applied to I and to a matrix B of constant vectors, the results are A" 1 and X where AX = B. Each 
transformation is selected so that the largest element is used in the pivotal position. (See ref. (a).) 

Accuracy: Total pivotal strategy is used to minimize the rounding errors; however, the accuracy of the 
final results depends upon how well- conditioned the original matrix is. 

Reference: (a) Fox, L.: An Introduction to Numerical Linear Algebra. Oxford Univ. Press, 1965. 

Storage : 542g locations. 
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A listing of subroutine MATIN V is as follows: 


SUBROUT IN E MATINVI A,N,B, M f OFTFKM, l PIVOT, INDEX, NMAX,T SCALE) 

********* OOf.lirtFNT DATE 08-01-6 8 SUBROUTINE REVISED 05-01-68 ********* 

MATRIX INVFftSIOM W I TH ACCOMPANYING SOLUTION OF LINEAR EQUATIONS 

0 I Nfc NS I (SM l P I V OT ( N ) , A ( NM AX » N ) ♦ BINMAX, M ) , l N OE X ( NMA X , 2 ) 

IQUI VALENCE ( IROlftt JROW), U C HI UM , J COUJM) t UMAX, T, SWAP) 

INI UAL 1/ AT I UN 

5 I$CALt=0 

6 R L r l 0, O** 10 J 

7 P.2- 1,0/P l 
10 nf TF«M*1.0 
15 D 1 20 J = L , N 
2 ) 1 PI VOT l J ) -0 
30 1)0 5 6 0 l = 1,\ 

S L a P CM r OK p) V 01 L;L FM ENT 

40 A 1 A X - U - 0 


4 6 0'1 

l<^ .1=1, 'i 




5 0 IF 

iii’iv'inji-ii 

4 0 » 

105 , 

30 

6 0 OJ 

UC K = 1 , N 




7 0 [F 

( r pi vnr< k ) - i i 

DO, 

1 00 , 

740 


8 0 IE I AHS( Afo/iX J-ABS ( M J,K ) ) ) 85 , 100 , 100 
8 5 IKC'WsJ 
90 1CL‘LU« = K 
9!) AMAX = A(J,<] 

1U0 CUKTIN'JF 
10 8 CUM'IM'jr 

IF UMAX) 110,106,110 
106 OE TERM =0,0 

isc.alf=o 
GO T 1 74 0 

110 l PI VOT ( 1C HUM) - I PIVOT 1 I C D L U M ) H 

c 

C IN TF PCrt /UT,r Pfho T P PUT PIVOT ELEMENT ON DIAGONAL 

C 

no IF UPLU-ICULUM) 140, 260, 140 

14 0 i;E T I f M= — UETFKM 

15 0 n l ?oi) (.■= 1 , N 

16 0 S V A P = A < I f ; : ! } ,v j L ) 

170 M IPJ-V ,UMI ICl:LUM,U 
200 A ( JCiJL MM ,|. ) =SWAP 
20 5 IFI M ) 2c0, 2 60, 210 
210 DO 2 60 l - 1 , r 

22 U SW AP = «‘. ( [*?*,{ ) 

23 0 0 U P'H.; , i. ) l C'.iLLJ H , I 1 
?60 f M I C ML iJM r L ) = SWAP 

260 I .JOE X (I , 1) = I PUW 

27 0 INDEX I 1,21= [ClJLUM 

5io pi vm = A( loauM, koi um» 

IE IPIVMT) luoo, 106, 1000 
r 

c S6/.U- THE DETERMINANT 

C 

luoo Pivni-iM/uT 

100 5 I >■ I A B S t ML T ? RM ) - rt l ) lu 30, 101 J, 1 0 L 0 
LOIG N ( “ r(-. * = ! u Tf-KM/o L 
l SCALES 1 $( Al t +• 1 

IF ( A K 5 ( I U TriR'M-Rl )lu60, 1020, 10 20 
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10? 0 L 

I SCALE- I SCALC+l 
Gil m ipaj 

1 03 0 I F ( A:^S( fH TERM I -R 2 ) 1040 , 1040 , 1060 
1 J 40 PlTLR v I-;)£T 

1 SC AL C = 1 SC M F-l 

!f (ArtSinirCEPM)-*?) 1 050 1 1 05 Ot 1060 
105 0 PfTt Kf' = nr TEPM*iU 
I SCALF= ISC ALC-1 

L 06 J IF (AMS(PIV)T IJ-RU 1090il070t 1070 

io7o pi vorr -pi vht i/m 

I SC.ALC= ISCALr+l 

!FIA'3S< PIVUT [ l-tU>320, 10 90 t 1 0 50 

i o*o p l vp it i =pi voti/ri 

I SCALP = I SCALOl 
r;n in 3?o 

1 090 I F ( AliM P I V-JT 1 I-R2 ) 2000 , 2000, 320 
2000 Pivnr I - P 1 Vi IT I*RL 
l S C A L F = I S C A L K — L 

I f ( Ar!S(P £ Vi it II-R2 120 10 i 2010, 3 20 
2010 PI VOTl=PI VHT !*RL 
I S (’ A L Fr = TSCALf - I 

3 ?o r>'UK^=DtrtRw#PivnTi 
C 

C D T V I OF PIVCT Ftp* iiV PIVOT ELEMENT 

C 

330 A( 1C0LUM, ICOL'JM >* I .0 
340 rrj 3 50 L = L * > J 

35 0 AUCOL'I^tL) = A{ ECOLUM,L I/P I VHT 

35 5 [MMJ 3fi0, 3 60, 360 

36 0 !) !”: 3 70 L-i ,m 

37o sm ir. ! .)nr-su = B( icoLUM t u/P! vnr 

PLHJCii i^niM-PI VUT ROWS 

3 90. in 550 ll=l,N 
390 IF ( 1. l-ir.ULJ«l 400, 550, 400 
400 T - A ( L 1 , rC’IUJM ) 

42 0 AlLl. [CTUJMl-0,0 

43 0 0 ) 4 60 L=l,N 

45 0 A I L l , I.) = A( L l *U-A ( ICliLU«t L ) *t 

45 5 I'-<N> r 50, 550 , 4 60 

46 0 P > 500 L = l , V 

60 0 rtlLl iL >“tHLl,L)-*( ICOLU^iL )*T 
5 60 CFMP'iir 

INTERCHANGE CCLlMNS 

600 DM 710 T -l ,\ 

610 L "= N ♦* 1 — T 

620 IF ( I Ni> t X ( L, D-INOtXlL t2M 630, 7L0, 

63 0 jrcwmn'h x a , 1 ) 

640 JCPL J M -IN0EX(L,2I 
65 0 r )n 705 K=1,N 
660 SwAP= A ( K, JROW) 

670 A(K, JROW) -A I K, JCOLUMI 
700 A <K,JC')LU-1 INSWAP 
706 CONTINUE 
710 CONTINUE 
740 RETURN 
END 
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STRAIN AND DISPLACEMENT TRANSDUCER LOCATIONS 
ON TORSION SPECIMEN 

Most torsion tests on tubular specimens are conducted with one or more strain 
rosettes and four LVDT's or similar instruments for measuring rotation. In order to 
keep user instructions at a minimum, a pattern for identifying the strain gages and 
rosettes was selected. (See fig. 10.) Starting with the first strain rosette and continuing 
clockwise around the circumference of the cylinder, the strain transducers are numbered 
in monotonically increasing order. The first and last transducers in each rosette are 
oriented at -45° and +45°, respectively, to the longitudinal axis of the specimen. As 
shown in figure 10, LVDT’s number 1 and 2 measure displacements of radial arms at a 
station near the rotated end of the tube (5^ and 62 ) whereas LVDT’s 3 and 4 measure 
displacements 63 and 64 , respectively, near the fixed end of the tube. Utilization of a 
different pattern for identifying transducers will generally require changes in the program 
in order to generate accurate and usable plots. 
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TABLE L- PROGRAM CAPABILITY 


Case 

Specimen 

Test 

Plot on Y -axis 

Plot on X-axis 

Compute slope 

Compute tangent modulus 

1 

Flat 

Compression 
or tension 

Q 

11 

%\* 

e l> € 2’ ’ ' *’ e 18 

No 

No 

2 

Flat 

Compression 
or tension 


e l> eg, . . anc * 

average of strains from up 
to four gages 

Yes, for averaged 
curve 

Yes, for averaged curve 

3 

Flat 

Compression 
or tension 


Average of strains from* up 
to four gages 

Yes 

Yes 

4 

Tube 

Compression 
or tension 

p 

ir(O.D. - t)t 

£ 1» f 2’ ’ • •’ e 18 

No 

No 

5 

Tube 

Compression 
or tension 

Q 

II 

o 

p* 

r+ 

~ — * 

r^l 

^2* ' * *j ^18 

average of strains from up 
to four gages 

Yes, for averaged 
curve 

Yes, for averaged curve 

6 

Tube 

Compression 
or tension 

+-> 

Z*' 

1 

6 

tl 

b 

Average of strains from up 
to four gages 

Yes 

Yes 

7 

Tube 

Torsion 

r 2T 

tt(O.D. - t) 2 t 

e l> e 2 f 18 

No 

No 

8 

Tube 

Torsion 

2T 

ir(O.D. - t) 2 t 

F 1 ,e 2 >’ * m € xB and sum of 
absolute values from ±45° 
gages in one rosette 

Yes, for sum of 
absolute values 
from ±45° 
gages 

Yes, for sum of absolute 
values from ±45° 
gages 

9 

Tube 

Torsion 

2T 

Average of sum of absolute 
values from all ±45° gage 

Yes 

Yes 

~ ir(O.D. - t) 2 t 

10 

L 

Tube 

Torsion 

2T 

Average of sum of absolute 
values from all ±45° gages 
and 

180(6! + 6 2 - 63 - 6 4 ) 

6 " 2 jtR£ 

Yes, for averaged 
strain curve 
and for torque - 
twist curve 

Yes, for averaged 
stress -strain curve 

ir(O.D. - t) 2 t 
and 

T 








Figure 2.- Program BECKPLT deck setup. 


58 











TANGENT MODULUS ,PSI 



h- + 


600 


400h + 


200 


oL 

o 


J L 


J L 


200 400 600 800 1000 1200 1400 1600 1800 

STRESS .PS I 


x 10 1 


Figure 5.- Tangent- modulus — compressive -stress curve developed from 
compression test on graphite -epoxy tube. 
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Figure 7.- Shear-stress— shear-strain curve 
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from torsion test on a graphite -epoxy tube 
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Overview of torsion specimen 
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Circumferential strain gaged section of specimen 


Figure 10,- Overview of torsion specimen and identification of strain-gage 

and LVDT numbering sequence. 
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